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Abstract 

This text is a slightly expanded version of my 6 hour mini-course at 
the PDE-meeting in Evian-les-Bains in June 2009. The first part gives 
some old and recent results on non-self-adjoint differential operators. 
The second part is devoted to recent results about Weyl distribution 
of eigenvalues of elliptic operators with small random perturbations. 

Resume 

Ce texte est une version legerement completee de mon cours de 
6 heures au colloque d'equations aux derivees partielles a Evian-les- 
Bains en juin 2009. Dans la premiere partie on expose quelques 
resultats anciens et recents sur les operateurs non-autoadjoints. La 
deuxieme partie est consacree aux resultats recents sur la distribution 
de Weyl des valeurs propres des operateurs elliptiques avec des petites 
perturbations aleatoires. 
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1 Introduction 

For self-adjoint and more generally normal operators on some complex Hilbert 
space "H we have a nice theory, including the spectral theorem and the well- 
known and important resolvent estimate, 

\\{z-Pr^<{d\st{z,a{P)))-\ (1.1) 

where cr(P) denotes the spectrum of P. The spectral theorem also gives very 
nice control over functions of self-adjoint operators, so for instance if P is 
self-adjoint with spectrum contained in the half interval [Aq, +oo[, then 

||g-tP|| < g-Aot^ t > 0. (1.2) 

However, non-normal operators appear frequently in different problems; 
Scattering poles. Convection-diffusion problems, Kramers- Fokker-Planck equa- 
tion, damped wave equations, linearized operators in fluid dynamics. Then 
typically, || {z — P)~^\\ may be very large even when z is far from the spectrum 
and this implies mathematical difficulties: 

- When studying the distribution of eigenvalues, 

- When studying functions of the operator, like e~^^ and its norm. 

The largeness of the norm of the resolvent far away from the spectrum 
also makes the eigenvalues very unstable under small perturbations of the 
operator and this is a source of mathematical and numerical difficulties. 

There are two natural reactions to this problem: 
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• Change the Hilbert space norm to make the operators look more nor- 
mal. This is quite natural to do when there is no clear unique choice 
of the ambient Hilbert space, like in problems for scattering poles (res- 
onances) . 

• Recognize that the region of the 2;-plane where — P)~^|| is large, has 
its own interest. One can then introduce the e-pseudospectrum which 
is the set of points 2; G C, which are either in the spectrum or such 
that 11(2; — -P)~"^|| > e. and study this set in its own right. 

The first point of view will be illustrated in the first part of these notes, 
and is at the basis of a whole range of methods from that of analytic ditations 
in the study of resonances to more microlocal methods. 

The second point of view has been promoted by numerical analysts like 
L.N. Trefethen and then made its way into analysis through contributions by 
E.B. Davies, M. Zworski and others. That spectral instability is not only a 
nuisance, but can be at the origin of nice and previously unexpected results, 
will hopefully be clear from the second and main part of these lecture notes, 
where we shall describe some results about Weyl asymptotics of the distri- 
bution of eigenvalues of elliptic operators with small random perturbations. 
These results have been obtained in recent works by M. Hager [30], [311 132] . 
Hager and the author |33], W. Bordeaux Montrieux |9], the author [881 ES] 
and Bordeaux Montriex and the author jlOj . 

The results in the second part of the lectures rely on microlocal analy- 
sis, combined with quite classical methods for non-self-adjoint operators and 
holomorphic functions of one variable and some elementary probability the- 
ory. It therefore was natural to include a first part that treats some older and 
newer results and methods about non-self-adjoint operators. Many of these 
methods combined with microlocal analysis have been used with great succes 
in resonance theory, but that is beyond the scope of these lectures.. Let us 
nevertheless mention results by T. Christiansen [HI [15] and Christiansen- 
P. Hislop [16] that establish Weyl type lower bounds for the number of scat- 
tering poles in large discs in generic situations. 

Here is the plan of the notes: 

Part I is devoted to some old and recent general results for non- self-adjoint 
differential operators. 

In Section [2] we describe some classical results for the distribution of 
eigenvalues of elliptic operators, starting with a result of T. Carleman about 
Weyl asymptotics of the eigenvalues for operators with real principal sym- 
bol. We also give a result by S. Agmon about the completeness of the set 
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of generalized eigenvalues as a well as further results by M.S. Agranovich, 
A.S. Markus, V.I. Matseev. 

In Section [31 we start by giving some basic definitions and facts about the 
e-pseudospectrum, then go on by describing the Davies-Hormander quasi- 
mode construction under the assumption that a certain Poisson bracket is 
non-zero and of suitable sign. 

In Section HJ we describe some estimates on the size of norm of the resol- 
vent when the spectral parameter is close to the range of the semi-classical 
principal symbol. These estimates are closely related to subellipticity esti- 
mates for operators of principal type and the results here are due to Dencker- 
Sjostrand-Zworski with some recent partial improvement of the author. 

In Section 0, we discuss some recent results from two different areas; 
the Kramers- Fokker- Planck operator and spectral asymptotics for analytic 
operators in two dimensions. The first topic (based on joint works with 
F. Herau-C. Stolk, F. Herau-M. Hitrik and inspired by works of Herau- 
F. Nier, B. Helffer-F. Nier) was chosen because it has a very concrete im- 
portance, while the second topic (based on joint works with M. Hitrik and 
S. Vii Ngoc) is of interest as an example of how to get precise information 
about individual eigenvalues. In both situations the methods exploit suitable 
changes of the Hilbert space norms. 

Part II is devoted to Weyl asymptotics for the eigenvalues of elliptic 
operators with small random perturbations. 

In Sectional we give a result, that generalizes and improves earlier results 
by Hager and Hager-Sjostrand about the number of zeros of holomorphic 
functions of exponential growth, which is clearly close to classical results for 
entire functions but that we have not found in the literature. This result is 
used in an essential way in the sequel. 

In Section [7] we treat the one- dimensional semi-classical case very much 
in the spirit of Hager. This case is easier than the general case, and it has 
some special features that permit to have more precise results, and is most 
likely the first testing case for more refined questions about statistics and 
correlation of eigenvalues. 

In Section Owe establish Weyl asymptotics in the multi-dimensional semi- 
classical case, by combination of complex analysis (Section E]), microlocal 
analysis spectral theory and probabilistic arguments. 

In Section [9|, we consider the large eigenvalues of elliptic operators. In 
the semi-classical case the results say that we have Weyl asymptotics with 
a probability tending to 1 very fast when Planck's constant tends to zero. 
The study of large eigenvalues of elliptic operators can often be reduced to 
a semi-classical study, and by performing such a reduction and applying the 
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Borel-Cantelli lemma, we show that Weyl asymptotics holds almost surely. 
Such a result was obtained by W. Bordeaux Montrieux for elliptic operators 
on using the results and the approach of Hager. Here we mainly describe 
a corresponding multi-dimensional result, obtained jointly with Bordeaux 
Montrieux, where the semi-classical part is the one described in Section [HI 
In Section (TD] we formulate some open problems. 



Part I 

Some general results 

2 Elliptic non-self-adjoint operators, some clas- 
sical results 

This is a classical area. T. Carleman [T2] considered the Dirichlet realization 
P of a second order elliptic operator in a bounded domain Q d R^, assum- 
ing enough smoothnes on the coefficients and on the boundary, and he also 
assumed that the principal symbol is strictly positive so that the non-self- 
adjointness can come only from the lower order symbols. In this case we see 
easily that the spectrum consists of isolated eigenvalues of finite algebraic 
multiplicity and we will allways count the eigenvalues with their multiplicity. 
Carleman showed that the eigenvalues are contained in a parabolic neighbor- 
hood of the positive real axis and that the real parts are distributed according 
to the Weyl asymptotics, i.e. 

e a{Py, sR/i < A} = i-A3/2(vol{(x,0 G T*n; < l}+o(l)), A ^ oo 

(2.1) 

His method of proof consisted in studying the trace of (P + k^) ^ — (P + Kq) ^ 
for a fixed kq in the limit when k — >■ oo and to apply a Tauberian argument. 

After Carleman there have been important results of Keldysch which 
have inspired later workers in the field, like Agmon, Agranovich, Markus 
and Matseev. In a spirit similar to that of Carleman, M.S. Agranovich and 
A.S. Markus [3] considered a non-self-adjoint elliptic classical pseudodiffer- 
ential operator on a compact manifold Q of dimension n, of order m > 0. 
Let p denote the principal symbol. Assume that the range of p is in a sector 
{z G C; |argp| < 6} where 9 < n. If the quantity 

'^'■=7^11 Pi^^O-'^^'SidOdx (2.2) 
(27r)'^n JnJ\(\=i 
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is 7^ then the authors show that 



N{X) X A"/™, A ^ oo 



(2.3) 



where iV(A) = #{(t{P) fl D{0, A)) and ^(O, A) denotes the open disc of ra- 
dius A centered at 0. They also obtain some more precise inequahties for 
the quantities hmsup;^^;.^^ A~"/™A^(A) and hminfA^>oo A~"/™'A^(A). Especially 
when these two limits are equal, then the corresponding quantity is in the 
interval A], where 



Here it is interesting to notice that d ^ when the angle 26 is smaller than 
nm/n. To prove these results, the authors first establish the asymptotic trace 
formula 



for im > n, where the constant can be expressed with the help of d, and then 
apply a Tauberian argument. Here {fi + A)~^ is of trace class when im > n. 

The operator P = f{x)Dx + g{x) on 5*^ is elliptic when f ^ 0, and 
its spectral behaviour can be studied explicitly. When the range of / is in 
a sector of angle less than vr, then we have a nice spectrum sitting on a 
line, while for larger values, strange things may happen ([H], [SS]) and the 
spectrum may be either equal to C or empty. 

Another interesting and potentially important question is that of the 
completeness of the set of generalized eigenfunctions. This question was 
studied by Keldysh and Agmon. Agmon studied elliptic boundary value 
problems, let us here formulate his result in the case of elliptic differential 
operators on manifolds without boundary [T] (see also [2]). Assume that 
P is such an operator of even order m > and assume that the symbol p 
takes its values away from a finite union of closed half-rays e*^^ [0, +oo[, where 
9i < 02 < ■■■ < On belong to [0, 27r[ and assume that the angles 9j+i — 9j are 
all strictly smaller than mir/n (with the convention that ^at+i = + 27r). 
Then the generalized eigenfunctions of P are complete in L^, in the sense 
that they span a dense subspace. 

The proof uses the following ingredients: 

• First we show that the resolvent {z — P)^^ is well-defined and of norm 
(9(|2;|^^) when z tends to infinity along one of the half rays given by 
aigz = 6j. In particular the resolvent exists for at least one value of z 
and by analytic Fredholm theory this implies that the spectrum of P 
consists of isolated eigenvalues and each such value is of finite algebraic 
multiphcity. 




(2.4) 



tr((/i + A)"^) ~ Const. yu' 



n/m—£ 



, n — )■ +00 



(2.5) 
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• Using a suitable functional determinant which is an entire function of 
fractional exponential grow and whose zeros are the eigenvalues, we 
get a polynomial control over the number A^(A) when A — oo and 
a corresponding exponential control over the resolvent on a family of 
circles dD{0,rj), rj — )■ oo, when j — >■ oo: 

<0(l)exp(|^|^+^), \z\=p,, (2.6) 

where e > can be chosen arbitrarily small. 

• Of course, {z — P)^^ will have poles at the eigenvalues, but if /o G 

is orthogonal to all the generalized eigenfunctions, then Ff{z) := 
{{z — -P)~^/|/o)l2 turns out to be an entire function for every / G L^. 
Thanks to the condition on the angles, and the exponential control of 
the resolvent, we can apply the Phragmen Lindelof theorem to conclude 
that Ff{z) is actually constant. By restricting to a ray of minimal 
growth we see that the constant has to be zero, and varying / it is not 
hard to see from this that /o has to be zero. 

Let us mention that rays of minimal growth of the resolvent have also 
been used by R. Seeley [82]. We could also point out that under the same 
ellipticity condition and the assumption that the range of p is not equal to 
all of C, we have the classical upper bound 

N{\) = C(A"/'"), A -> +00. (2.7) 

This can be proved in a similar manner using a relative determinant and 
the Jensen formula. A more elementary proof is the following, that was 
pointed out to me by M.S. Agranovich [1]: After replacing P by P — Aq for a 
sufficiently large Aq on a ray of minimal growth, we may assume that P is bi- 
jective. Let si(P"^) > S2(P"^) > ... be the singular values of the inverse P 
(i.e. the decreasing sequence of eigenvalues of (P*P)~^/^ which is bounded 
from L2 to H""). If C > is large enough we have {P*P)-^ < C{1 - A)""* 
where A is the Laplace Beltrami operator on X for some smooth Riemannian 
metric and this implies corresponding inequalities for the eigenvalues. Ap- 
plying the Weyl asymptotics for —A, we conclude that Sj{P~^) < 0{j~"^^"'). 
Then Corollary 3.2 in Chapter III of [27] implies that A^fp^^) = C(j-™/"), 
where \j{P~^) denote the eigenvalues of P~^ arranged so that j i— )■ |Aj(P~^)| 
is decreasing. Then (12. 7p follows. 

Markus and V.I. Matseev [53] have established interesting estimates on 
the difference of the counting function for a self-adjoint (or more generally 
normal) operator and the counting function for the real parts of the eigen- 
values of a small perturbation of that operator. The proofs are based on 
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the use of relative determinants. We will not state the general results, but 
simply mention a corollary of the general result: Let P be an elliptic dif- 
ferential operator of order m on a compact manifold with positive principal 
symbol ^). Then the eigenvalues are contained in a thin neighborhood 
of the positive real axis and if A^(A) denotes the number of such eigenvalues 
with real part < A, then A^(A) = (27r)^"Volp^i([0, 1])A"/"' + C(A("-i)/™). 
Notice that the remainder estimate is the same as in the general result of 
Avakumovic [6], Levitan [62] (m = 2) and Hormander [39] (general m) for 
self-adjoint elliptic operators which in turn depend on trace formulas, not for 
resolvents as in Carleman's approach but for hyperbolic evolution problems. 

3 Pseudospectrum, quasi-modes and spectral 
instability 

Let "H be a complex Hilbert space and let P : "H — "H be a closed densely 
defined operator. Recall that the resolvent set is defined as 

p{P) = {-2 G C; P — 2; : T>{P) "H has a bounded 2-sided inverse}. 

It is an open set, if 2; G p(P) and \\{z — P)^^|| = 1/e, then the open disc 
D{z, e) is contained in p(P). The spectrum of P is the closed set 

a(P) = C\p(P). 

Following Trefethen-M. Embree [HZ] we define, for e > 0, the e-pseudospectrum 
as the open set 

a,(P) = a{P) U{ze p(P); - P)-'\\ > 1/e}. (3.1) 

Unlike the spectrum, the e-pseudospectrum will change if we replace the 
given norm on H by an equivalent one. 

<Je{P) can be characterized as a set of spectral instability, by the following 
simplified version of a theorem of Roch and Silberman: 

a,(P)= U a(P + Q). (3.2) 

IIQII<e 

The result becomes more sublte if we use the more traditional definition with 
a non-strict inequality in f l3.ip . 

Proof. Let 5e(P) denote the right hand side in (13. 2p . If z G C \ a^{P), then 
by a perturbation argument, we see that z G C \ o^e{P)- 
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Let z G aSP\ If 2; G cr(P) we also have z G (^^{P), so we may assume that 
z G p{P). Then 3 n G G "H such that ||n|| = 1, < e, {P-z)u = v. 

Let Q be the rank one operator from "H to "H, given by (50 = — (0|m)v. Then 
= ll'^llll'^ll < e and (P + Q-2)m = v + Qu = v-v = 0, so z e a{P + Q), 
and 2; G ae(-P)- D 

Using the subharmonicity of the function z t— ?■ — P)~^|| we notice that 
every bounded connected component of (Je{P) contains an element of cr(P). 

We next discuss the construction of quasimodes for non-normal differ- 
ential operators which shows that very often we get large e-pseudospectra. 
The background and starting point is a result by E.B. Davies [19] for non- 
selfadjoint Schrodinger operators in dimension 1. M. Zworski |103] observed 
that this is essentially an old result of Hormander |171 |1H](1960), and that 
we have the following generalization, with {a, h} = a'^ ■ h'^ — a'^ ■ h'^ = Haip) 
denoting the Poisson bracket of a = a(x,^), 

Theorem 3.1 Let 

P{x,hD,) = J2 aa{x){hD,T, = (3.3) 

\a\<m 

have smooth coefficients in the open set Vt C R". Putp{x,^) = J2\a\<m. (^oi{x)^' 
Assume z =p{xo,^o) with ]{p,p}{xo,^o) > 0. Then 3 u = Uh, with \\u\\ = 1, 
II (P - z)ti\\ = 0{h°^), when h^O. 

In the case when the coefficients are all analytic we can replace by 

U^-l/Ch gQJ^g C > 

Notice that this implies that if the resolvent {P — z)~^ exists then its norm 
is greater than any negative power of h when /i — > (and even exponentially 
large in the analytic case). 

In the case n > 2, we noticed with A. Melin in [70] that ii z = p{p) 
and ^p, '^p are independent at p, then j{p,p} times the natural Liouville 
measure is equal to a constant times the restriction to p~^{z) of o"""^ which 
is a closed form. It follows that if F is a compact connected component of 
p~^{z) on which d^p and d'^p are pointwise independent, then the average of 
over r with respect to the Liouville measure has to vanish. Hence if 
there is a point on F where the Poisson bracket is 7^ then there is also point 
where it is positive. In the case n = 1 we have a similar phenomenon: If 
(for instance thanks to suitable ellipticity assumption) we know that p~^{z) 
is finite and that \{p,p} is 7^ everywhere on that set, then this set is finite 
and if it is contained in the interior of a connected bounded set Q in phase 
space with smooth boundary such that the variation of arg {p — z) along that 



10 



boundary is equal to zero, then we have to have an equal number of points in 
p~^{z) where ^{p,p} is positive and where it is negative. This follows from 
the observation that the argument variation of p — z along a small positively 
oriented circle around a point in p~^{z) is =p27r when ±^{p,p} > at that 
point. 

Example 3.2 P = -h'^A + V{x), p{x,C) = + V{x), ]{p,p} = -4^ ■ 
QV'{x). 

More recently K. Pravda-Starov [71] improved this result by adapting 
a more refined quasi-mode construction of R. Moyer (in 2 dimensions) and 
Hormander [51] for adjoints of operators that do not satisfy the Nirenberg- 
Treves condition for local solvability. 

The proof in the C°°-case in |1U3] is by a standard reduction of semi- 
classical results to classical results in ordinary microlocal analysis. In [22] we 
gave a direct proof and also treated the case of analytic coefficients, which is 
also essentially quite old. Here is a brief outline of a 

Proof. In the following we use the notation neigh (a. A) for "some neigh- 
borhood of a in A\ If G C~ (neigh (xq, R")) satisfies (j)'{o) : e R", 
and 

Q(P"ixo) > 0, (3.4) 
then we can define the complex Lagrangian manifold 

A,^ := {(x, (P'ix)); x € neigh (xq, C")} (3.5) 

where we extend to a complex neighborhood by taking an almost holomor- 
phic extension, i.e. a smooth extension such that (90 = 0{{'^x)°°). In this 
case we can content ourselves with working with formal Taylor expansions at 
Xo, and then can be viewed as an equivalence class of real submanifolds 
of the complexified phase space C^" where two submanifolds are equivalent 
if they agree to infinite order at {xq, $,q). As observed by Hormander [50] and 
developped a lot by the author with A. Melin in [BH] and in other works, the 
positivity assumption 03.41) can be formulated equivalently by saying that 

^a{t,t)>0, OT^t gT(,„,5„)(A^), (3.6) 

where a denotes the symplectic 2-form, here viewed as a bilinear form on the 
complexified tangentspace of the cotangent space at (xo,^o)- 

Let z, p, (xo, ^o) be as in the theorem. Then we observe that \a{Hp, Hp) = 
^{p,p} > 0. Moreover, the real set S := p~^{z) is a smooth symplec- 
tic manifold near (xq, i^o) and using the Darboux theorem, we can identify 
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it with R2("~i) and hence find a Lagrangian submanifold A' in its com- 
pexification passing through (xc^o) that satisfies the positivity condition 
(13.61) . Viewing the complexification of S as a submanifold of C^", we can 
take A = {exp sHp{p); s G neigh (0,C),p G neigh ((xq, ^o), C^")}. Using 
that Hp is symplectically orthogonal to the tangent space of S it is then 
quite easy to verify that A is a complex Lagrangian manifold to oo order 
at (xoi^o) contained (to infinite order) in the complex characteristic hyper- 
surface {p G neigh ((xo, ^o), C^""); p(p) = 0} and satisfying the positivity 
condition (13.61) . Hence to infinite order, A is of the form for a function (p 
as in (13. 4p . (13. 5p . which also fulfills the eiconal equation 

p(x,0'(x)) = O(|x-Xor). (3.7) 

We normalize by requiring that 0(xo) = 0. Then the function e^'t'i^)/'^ 
is rapidly decreasing with all its derivatives away from any neighborhood 
of Xq, and by a complex version of the standard WKB-construction we 
can construct an elliptic symbol a{x;h) x ao(x) + hai{x) + by solv- 
ing the suitable transport equations to infinite order at Xq, such that if x ^ 

(neigh (xq, R-")) is equal to 1 near xq, then u(x; h) = x{x)h~"'^^a{x; h)e^'^^^^^'^ 
has the required properties. □ 



4 Boundary estimates of the resolvent 
4.1 Introduction 

In this section we are interested in bounds on the resolvent of an /i-pseudo- 
differential operator when z is close to the boundary of the range of p. As with 
the quasi-mode construction this question is closely related to classical results 
in the general theory of linear PDE, and with N. Dencker and Zworski (|22]) 
we were able to find quite general results closely related to the classical topic 
of subellipticity for pseudodifferential operators of principal type, studied by 
Egorov, Hormander and others. See [5T] . 

In [22] we obtained resolvent estimates at certain boundary points, 

(A) under a non-trapping condition, 
and 

(B) under a stronger "subellipticity condition" . 

In case (A) we could apply quite general and simple arguments related 
to the propagation of regularity and in case (B) we were able to adapt gen- 
eral Weyl-Hormander calculus and Hormander's treatment of subellipticity 
for operators of principal type ([51]). In the first case we obtained that 
the resolvent extends and has temperate growth in 1/h in discs of radius 
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0{h\nl/h) centered at the appropriate boundary points, while in case (B) 
we got the corresponding extension up to distance 0{h''^^''~^^^), where the 
integer k > 2 is determined by a condition of "subelhpticity type" . 

Using a method based on semi-groups led to a strengthened result in case 
(B): The resolvent can be extended to a disc of radius 0({h\nl/h)''^^''^^^) 
around the appropriate boundary points. 

Let X be equal to R" or equal to a compact smooth manifold of dimension 

n. 

In the first case, let m G C°°(R^"; [1, +oo[) be an order function (see [23] 
for more details about the pseudodifferential calculus) in the sense that for 
some Co, Nq > 0, 

m{p)<Co{p-fif''m{fi), p,/ieR2", (4.1) 

where (p — //) = (1 + |p — pp)^/^. Let P = P(x,^; h) E S{m), meaning that 
P is smooth in x, and satisfies 

\d^^^P{x,^;h)\ < C„m(x,0, (x,0 e R2", a e N^", (4.2) 

where Ca is independent of h. We also assume that 

Pix,^;h) po{x,0 + hpiix,0 + m^(m), (4.3) 

and write p = Po for the principal symbol. We impose the ellipticity assump- 
tion 

3w e C, C > 0, such that \p{p) - w\ > m{p)/C, Vp G R^", |p| > C. (4.4) 
In this case we let 

P = P'^ix, hD^; h) = Op(P(x, h^; h)) (4.5) 

be the Weyl quantization of the symbol P{x, hC,', h) that we can view as a 
closed unbounded operator on L^(R"'). 

In the second case when X is compact manifold, we let P G S'"o(T*X) 
(the classical Hormander symbol space ) of order m > 0, meaning that 

\d:dlP{x,i-h)\ < C^Air-^^\ {^.i) e T*X, (4.6) 

where Ca^p are independent of h. We also assume that we have an expansion 
of the type fl4.3p . now in the sense that 

Af-l 

P{x,i;h) -J^h'pA^^O e h^S^.-'^iT-'X), N = 1,2,... (4.7) 
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and we quantize the symbol P{x, hC,; h) in the standard (non-unique) way, 
by doing it for various local coordinates and paste the quantizations together 
by means of a partition of unity. In the case m > we impose the ellipticity 
condition 

3C > 0, such that > ^7^, \^\ > C. (4.8) 



Let S(p) = p*{T*X) and let Soo(p) be the set of accumulation points 
of p{pj) for all sequences pj G T*X, j = 1,2,3, .. that tend to infinity. By 
pseudo differential calculus, if K C C\S(p) is compact, then {z — P)~^ exists 
and has uniformly bounded operator norm for z & K, < h <^ 1. 

The following theorem ([90]) is a partial improvement of corresponding 
results in 1221. 



Theorem 4.1 We adopt the general assumptions above. Let zq G 9S(p) \ 
Soo(p) (in-d assume that dp ^ at every point of p^^{zq). Then for every 
such point p there exists 6 & H ( unique up to a multiple of n) such that 
d{e~^^{p — 2o)) i'S real at p. We write 9 = 9{p). Consider the following two 
cases: 

• (A) For every p G p~^{zq), the maximal integral curve of H^^^-ie{p)p^ 
through the point p is not contained in p~^{zo). 

• (B) There exists an integer k > 1 such that for every p G p~^{zo), there 
exists j G {1, 2, .., /c} such that 

p*{exptHp{p)) = at^ + C(t^+^), t -> 0, 

where a = a{p) 7^ 0. Here p also denotes an almost holomorphic exten- 
sion to a complex neighborhood of p and we put p* {p) = p(jl) . Equiva- 
lently, Hl{p){p)/{j\) = a^O. 

Then, in case (A), there exists a constant Cq > such that for every 
constant Ci > there is a constant C2 > such that the resolvent {z — P)^^ 
is well-defined for \z — zq\ < Ci/iln^, h < and satisfies the estimate 

||(z-P)-i||<^exp(^|z-Zo|). (4.9) 

In case (B), there exists a constant Cq > such that for every constant 
Ci > there is a constant C2 > such that the resolvent {z — P)~^ is 
well-defined for \z — zq\ < Ci{h\\ij^Y/^^^^\ h < and satisfies the estimate 

||(.-P)-||<^exp(§^|.-.o|^). (4.10) 
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In [22] we obtained fl4.1Up for z = zq, implying that the resolvent exists 
and satisfies the same bound for \z — zo\ < h^^'^^^^'^ / 0{1) in case (B) and 
with k/{k + 1) replaced by 1 in case (A). In case (A) we also showed that 
the resolvent exists with norm bounded by a negative power of h in any 
disc D{zQ,Cihhi{l/h)). (The condition in case (B) was formulated a little 
differently in [22], but the two conditions lead to the same microlocal models 
and hence they are equivalent.) The case (A) of the theorem is basically 
identical with the corresponding result in [22] and was proved using weighted 
estimates with weights that have at most polynomial growth in h. We will 
not disuss that in detail here and instead we concentrate on the partially new 
result in the case (B). 

When k = 2 more direct methods are available and more precise bounds 
can be given, at least in special cases. Such results have been obtained by 
J. Martinet [66], Y. Almog, B. Helffer, X. Pan, see [31] and W. Bordeaux 
Montrieux [9]. 

Let us now consider the special situation of potential interest for evolution 
equations, namely the case when 

zo e iR, (4.11) 
3fJp(p) > in neigh {p-\zo),T*X). (4.12) 



Theorem 4.2 We adopt the general assumptions above. Let zq G dT,{p) \ 
Soo(p) o-'iT'd assume (^TTI]), (^T7^. Also assume that dp on p~^{zq), so 



that d'^p 7^ 0, d^p = on that set. Consider the two cases of Theorem\4.1 



• (A) For every p E p ^{zq), the maximal integral curve of Hc^p through 
the point p contains a point where ^p > 0. 

• (B) There exists an integer k > 1 such that for every p G p~^{zo), we 
have HiiJR,p{p) ^ for some j G {1,2, k} . 

Then, in case (A), there exists a constant Cq > such that for every 
constant Ci > there is a constant C2 > such that the resolvent {z — P)~^ 
is well-defined for 

1-1 11 
\'^{z-zo)\ <—, —<^z< Cih\n-, h<—. 



and satisfies the estimate 

^exp{^^z),^z > -h. 
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In case (B), there exists a constant Co > such that for every constant 
Ci > there is a constant C2 > such that the resolvent [z — P)~^ is 
well-defined for 

\^{z-z,)\<^, ^<3ft^<Ci(/ilni)^, h<^, (4.14) 
and satisfies the estimate 

4.2 Outline of the proof in case (B) 

Away from the set p~^{zo) we can use ellipticity, so the problem is to obtain 
microlocal estimates near a point po ^ P'^i^o)- After a standard factorization 
of P — z in such a region, we can further reduce the proof of the first theorem 
to that of the second one. 

The main (quite standard) idea of the proof of Theorem 14.21 is to study 
exp{—tP/h) (microlocally) for < t <^ 1 and to show that in this case 

tP t^'^^ 
llexp- — II <Cexp(— ^), (4.16) 

for some constant C > 0. Noting that that imphes that || exp — ^|| = 0{h°°) 
for t > when 5{k + 1) < 1, and using the formula 

{z-P)-' = --J^ exp{^-^)dt, (4.17) 

we get to fHTT^ . 

The most direct way of studying exp{—tP/h), or rather a microlocal 
version of that operator, is to view it as a Fourier integral operator with 
complex phase ([671 EHl EHl [68]) of the form 



^^^^^"^^ = j j a;, rj; h)u{y)dydr,, (4.18) 

where the phase should have a non-negative imaginary part and satisfy the 
Hamilton- Jacobi equation: 

idt^ + p{x, d^cfy = C((S>0)°°), locally uniformly, (4.19) 



16 



with the initial condition 



0(0, X, 77) = X ■ 77. (4.20) 

The amphtude a will be bounded with all its derivatives and has an asymp- 
totic expansion where the terms are determined by transport equations. This 
can indeed be carried out in a classical manner for instance by adapting the 
method of [69] to the case of non- homogeneous symbols following a reduction 
used in [721 [68]. It is based on making estimates on the fonction 

S^it) = '^{[ Cis) ■ dx{s)) - ^^{t) ■ ^>x(t) + sRe(O) ■ '^x{0) 
Jo 

along the complex integral curves 7 : [0, T] 9 s i-t- (x(s),^(s)) of the Hamil- 
ton field of p. Notice that here and already in fl4.19p . we need to take an 
almost holomorphic extension of p. Using the property (B) one can show 
that '^(j){t,x,T]) > C~^t*^"*"^ and from that we can obtain (a microlocalized 
version of) f l4.16p quite easily. 

Finally, we prefered a variant: Let 



3n 



Tu{x) = Ch-- J e^^^''^y>u{y)dy, 

be an FBI - or (generalized) Bargmann-Segal transform that we treat in the 
spirit of Fourier integral operators with complex phase as in [85]. Here (j) 
is holomorphic in a neighborhood of (xo,yo) ^ C" x R", and — 0^(a;o,|/o) = 
r/o G R", Q(l)ly{xo,yo) > 0, det^'^^y{xo,yo) ^ 0. Let Kt : {y, -cj)'y{x,y)) ^ 
{x , (f)'^{x , y)) be the associated canonical transformation. Then microlocally, 
T is bounded H^^ := Hoi (fi) n L^{n, e~'^'^''/^L{dx)) and has (microlo- 

cally) a bounded inverse, where is a small complex neighborhood of xq in 
C". Here the weight $0 is smooth and strictly pluri-subharmonic. If : = 
{(x, 1^^); X G neigh (xo)}, then (in the sense of germs) = Kt{T*X). 
The conjugated operator P = TPT~^ can be defined locally modulo 0{h°°) 
(see also ^7\) as a bounded operator from — )■ provided that the 
weight $ is smooth and satisfies $' — $q = 0{h^) for some 5 > 0. (In the 
analytic frame work this condition can be relaxed.) Egorov's theorem applies 
in this situation, so the leading symbol p of P is given by p o = p. Thus 
(under the assumptions of Theorem 14. 2p we have 3f?P| > 0, which in turn 

can be used to see that for < t < , we have e~^^/^ = 0{1): H^^ — )■ H^^, 
where $t < $0 is determined by the real Hamilton- Jacobi problem 

d^t 2 d^t 
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Here is a somewhat formal derivation of ( I4.2ip : 
Consider formally: 



and try to choose $t so that the time derivative of this expression vanishes 
to leading order. We get 



« hdt j UtUte-'^^'/^L{dx) 

= - (^{Pu,\ut)H,^ + {ut\Put)H,^ + j 2^{x)\u\'e-^''^'^L{dx)^ . 

Here 

{Pnt\ut)H,^ = I {p\,^^+0{h))\ut\'e-'''^/'L{dx), 
and similarly for {ut\Put)H^^, so we would like to have 

^ I (2^ + 2^p\,^^ + Omiu.l'e-'-'^/'Lidx). 

We choose $t to be the solution of (14.21 1) . Then the preceding discussion 
again shows that e~^^/^ = (9(1) : if^^ H^^. 

To get dUSD, it suffices to show that < $0 - f'+^/C for < t < 1. 
By a geometric discussion, this follows from 

Gt{p) < -t'+VC (4.22) 
where Gt is a smooth function in a real neighborhood of po, given by 



dt 



^p{p + iHgAp)) = ^. ^0 = 0. (4.23) 



The behaviour of Gt is easy to understand by means of Taylor expansion of 
the first equation in 04.231) at the point p. 



4.3 Examples 

Consider 

P = -h'^A + iV{x), V E C°°(X; R), (4.24) 

where either X is a smooth compact manifold of dimension n or X = R". 
In the second case we assume that p = + iV{x) belongs to a symbol 
space S{m) where m > 1 is an order function. It is easy to give quite general 
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sufficient condition for this to happen, let us just mention that if V G C^(R^) 
then we can take m = 1 + ^'^ and if d°'V{x) = 0{{1 + for all a G N*^ 

and satisfies the ellipticity condition |V^(a;)| > C~^|xp for |x| > C, for some 
constant C > 0, then we can take m = 1 + + x"^. 

We have S(j9) = [0, oo[+iV {X) . When X is compact then Soo(p) is empty 
and when X = R", we have Soo(p) = [0, oo[+iSoo(^), where T^ooiV) is the 
set of accumulation points at infinity of V. 

Let zo = xo + iyo G \ Soo(p). 

• In the case Xq = we see that Theorem 14. 2 1 (B) is applicable with k = 2, 
provided that is not a critical value of V. This is close to problems 
from fluid dynamics, studied by I. Gallagher, T. Gallay, F. Nier [26] . 

• Now assume that Xq > and that is either the maximum or the min- 
imum of V. In both cases, assume that V^^{yo) is finite and that each 
element of that set is a non-degenerate maximum or minimum. Then 
Theorem 14.21 (B) is applicable to ±iP. By allowing a more complicated 
behaviour of V near its extreme points, we can produce examples where 

(B) applies with k > 2. 

Now, consider the non-self-adjoint harmonic oscillator 

Q = -^,+w' (4.25) 

on the real line, studied by Boulton [TTj and Davies [20], K. Pravda Starov 
[75] . Consider a large spectral parameter E = iX + fi where A ^ 1 and <^ 
A. The change of variables y = y/Xx permits us to identify Q with Q = XP, 
where P = -/i^^ + ix^ and /i = 1/A ^ 0. Hence Q-E = A(P-(i + f)) and 
Theorem 14.21 (B) is applicable with k = 2. We conclude that {Q — E)~^ is 
well-defined and of polynomial growth in A (which can be specified further) 
when 

Y <C'i(A^MnA)i, 
A 

for any fixed Ci > 0, i.e. when 

/i < CiA^(lnA)i (4.26) 

M. Hitrik and K. Pravda Starov [H] have obtained interesting results on 
the characterization of the exponential decay of the the semi-groups gener- 
ated by differential operators with quadratic symbols. 
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5 Survey of some recent results 



5.1 Introduction 

In this section we survey some recent results about Kramers-Fokker-Planck 
type operators and about Bohr- Sommerf eld quantization conditions in di- 
mension 2. In both cases our approach uses quite essentially the possibility 
of modifying the Hilbert space structure by means of exponential weights. 

In recent years there has been many works that apply the commutator 
methods developed by Kohn for subelliptic operators to equations of Fokker- 
Planck type and non-equilibriaum stat physics models. Especially there is 
the work by F. Herau-F. Nier [JT] devoted to the Kramers-Fokker-Planck 
operator 

P = yhd,-V\x)-hdy+'^{y-hdy)-{y+hdy), gR'" = R^xR;;. (5.1) 

Symbol: p = ip2 + Pii where 

^, ri) = \{y'^ + V^), P2{x, y,^,v) = y-^- V'{x) ■ r] (5.2) 

Here we have suppressed some physical parameters and only kept h which is 
proportional to the temperature. Using commutator techniques Herau and 
Nier establish an interesting link to a Witten Laplacian and under various 
general symbol type assumtions on V they show: 

• global subellipticity (P is not elliptic even locally), and m-accretivity, 

• absence of large eigenvalues and corresponding power-decay of the re- 
solvent in certain "parabolic" neighborhoods of zR, 

• Estimates on the smallest non-vanishing eigenvalue (relating it to the 
corresponding quantity for the Witten laplacian), and especially esti- 
mates on this quantity in the high and low temperature limit. 

• Precise estimates on the return to equilibrium (when is an eigenvalue) 
or simply decay when is not an eigenvalue, including estimates on the 
exponential rate of convergence. 

With F. Herau and C. Stolk |32] we applied microlocal methods in order to 
study the semi-classical (low temperature) limit, the results are fewer than 
the ones in [H] and for me easier to describe. 
Assume that ^ is a Morse function, such that 

|y'(a;)| > 1/C when > C, ^"^"(a;) = 0(l),Va G N". (5.3) 
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We showed that 3c > such that for any C > 1 we have have for /i > 
small enough: 

• The eigenvalues in the disc D{0,Ch) are of the form /i = /i(A) = 
Xh + o{h), where A are the eigenvalues of the quadratic approximation 
of P\f^_^ at the critical points of V. 

• For \z\ > Ch and < c|z|^/'^/i^/^ the resolvent exists and satisfies the 
estimate 

This is similar to Theorem 14.21 (B) in the case k = 2. with an essential 
difference: p has no values in iH \ {0}, so (9S C Sqo! Nevertheless we can 
compute Hp^pi = V'{xY + — (y"{x)y ■ y + Vi^x)!] ■ t)) and see that when 
Pi is small, then H^^pi > except near points where ^ = V'{x) = 0. This 
means that where pi is very small, the short time averages are not small. 
This can be exploited by conjugating the operator by an operator which is 
bounded and has a bounded inverse and for which the conjugated operator 
has a symbol with a larger real part. 

Below we shall discuss a more general supersymmetric situation and also 
give detailed results on the return to equilibrium for the associated heat 
equation, following two joint works with F. Herau and M. Hitrik. 

In the case of the KFP equation we make a moderate change of the Hilbert 
space norm in order to increase the real part of the operator away from certain 
critical points of the symbol. When we have analyticity assumptions such 
changes of the norm can be much larger and sometime allow us to study 
all eigenvalues also at fixed distance inside the range of p. Here is one such 
result by M. Hitrik |13]: 

Let P = P{x, hD^] h) on R satisfying the general conditions of section |H 
Also assume that P has a holomorphic extension to a tubular neighborhood 
of in which is still (9(m(3ft(x, ^)) Let zq G 9S \ Eqo be a point such 
that p~^{zo) = {{xq, ^o)} and such that \p{x,C,) — zo\ x |(x — xq, ^ — ^o)P for 
(x,^) in a neighborhood of {xo,^q). Also assume that there is a truncated 
sector -20+] 0, eo]e*[^"~'^"'^''"'"'^°l which is disjoint from Then in a small but 

fixed neighborhood of zq, the spectrum is given by the set of values zq + 
G{{k + h) + where G{-; h) is holomorphic and ~ Y.'^ h^Gj{-; h), 

k = 0,l, 2, where Go(0) = 0, G'^iO) ^ 0. 

The main idea is to construct an IR-manifold A (a complex deformation 
of real phase space), containing {xo,C,o), such that on A, p — zq is elliptic 
outside an arbitrarily small neighborhood of (xo,^o) and such that near that 
point, the restriction of the quadratic part of to A takes its values along a 
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ray in the complex plane. To implement this picture, we make a Bargmann 
transform, mapping into an exponentially weighted space of holomorphic 
functions, then deform the weight. 

In the second part of this section we discuss precise Bohr- Sommerf eld 
rules in dimension 2 for non-self-adjoint operators in the semi-classical limit. 
What is remarkable here is that thanks to the non-self-adjointness we get 
better results than what would be possible in the self-adjoint case. Again 
analyticity assumptions and the use of exponential weights on the Bargmann 
transform side is essential. We will follow recent works with M. Hitrik and 
S. Vn Ngoc gSl Sg 

5.2 Kramers-Fokker-Planck type operators, spectrum 
and return to equilibrium 

5.2.1 Introduction 

There has been a renewed interest in the problem of "return to equilibrium" 
for various 2nd order operators. One example is the Kramers-Fokker-Planck 
operator: 

P = yhd^- V'{x) ■ hdy + ^{-hdy + y) ■ {hdy + y), (5.1) 

where x, y G R" correspond respectively to position and speed of the particles 
and > corresponds to temperature. The constant 7 > is the friction. 
(Since we will only discuss aspects we here present right away an adapted 
version of the operator, obtained after conjugation by a Maxwellian factor.) 
The associated evolution equation is: 

{hdt + P)u{t,x,y) = Q. 

Problem of return to equilibrium: Study the rate of convergence of u{t,x,y) 
to a multiple of the "ground state" uo{x, y) = e~^^ /2+vix))/h -^j^g^ f _^ -|-oo, 
assuming that V{x) — t- +00 sufficiently fast when x — )■ 00 so that Uq G 
L^(R^"). Notice here that P{uq) = and that the vector field part of P is 
h times the Hamilton field of y'^/2 + V{x), when we identify R^"^^ with the 
cotangent space of R". 

A closely related problem is to study the difference between the first 
eigenvalue (0) and the next one, /i(/i). (Since our operator is non-self-adjoint, 
this is only a very approximate formulation however.) 

Some contributions: L. Desvillettes-C. Villani [21], J. P. Eckmann-M. Hairer 
[25], F. Herau-F. Nier [H], B. Helffer-F. Nier [35], Villani [98]. In the work 
[H] precise estimates on the exponential rates of return to equilibrium were 
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obtained with methods close to those used in hypoelhpticity studies and this 
work was our starting point. With Herau and C.Stolk [12] we made a study in 
the semi-classical limit and studied small eigenvalues modulo 0{h°°). More 
recently with Herau and M. Hitrik [39] we have made a precise study of 
the exponential decay of when V has two local minima (and in that 
case turns out to be real). This involves tunneling, i.e. the study of 
the exponential decay of eigenfunctions. As an application we have a pre- 
cise result on the return to equilibrium [ID]. This has many similarities with 
older work on the tunnel effect for Schrodinger operators in the semi-classical 
limit by B. Helffer-Sjostrand [361137] and B. Simon [HI] but for the Kramers- 
Fokker-Planck operator the problem is richer and more difficult since P is 
neither elliptic nor self-adjoint. We have used a supersymmetry observation 
of J.M. Bismut [7] and J. Tailleur-S. Tanase-Nicola-J. Kurchan [HS], allowing 
arguments similar to those for the standard Witten complex [37] . 

5.2.2 Statement of the main results 

Let P be given by ([5lD where V G C°°(R"; R), and 

a°V^(x) = Oil), \a\ > 2, (5.2) 

|VV^(x)| > 1/C, \x\ > C, (5.3) 

y is a Morse function. (5-4) 

We also let P denote the graph closure of P from iS(R^"^) which coincides with 
the maximal extension of P in (see [HI [351 13H])- We have 3fJP > and 
the spectrum of P is contained in the right half plane. In [12] the spectrum 
in any strip < < Ch (and actually in a larger parabolic neighborhood 
of the imaginary axis, in the spirit of [H]) was determined asymptotically 
mod {0{h°°)). It is discrete and contained in a sector \Qz\ < C^z + 0{h°°): 

Theorem 5.1 The eigenvalues in the strip < < Ch are of the form 

\,4h) ~ hifi.^k + h^'^^'-'^HKi + ^'^'^^'>j>,2 + ••) (5.5) 

where fij^k O'l"^ the eigenvalues of the quadratic approximation ( "non- self adjoint 
oscillator") 

2/ ■ <9x - V"{xj)x ■ dy + ^{-dy + y) ■ {dy + y), 
at the points {xj,0), where Xj are the critical points of V . 
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The known explicitly and it follows that when Xj is not a local 

minimum, then 5RAj fc > h/C for some C > 0. When Xj is a local minimum, 
then precisely one of the Xj^k is 0{h°°) while the others have real part > h/C . 
Furthermore, when V — ?■ +00 as x — )■ 00, then is a simple eigenvalue. In 

particular, if V has only one local minimum, then 

inf K(a(P) \ {0}) ~ + h^^2 + ...), /ii > 0. 

(or possibly an expansion in fractional powers) and we obtained a corre- 
sponding result for the problem of return to equilibrium. It should be added 
that when simple eigenvalue of the quadratic approximation then 

Nj^k = 1 so there are no fractional powers of h in (15. 5p . 

The following is the main new result that we obtained with F. Herau and 
M. Hitrik in |39]: 

Theorem 5.2 Assume that V has precisely 3 critical points; 2 local minima, 
x±i and one "saddle point", xq of index 1. Then for C > sufficiently large 
and h sufficiently small, P has precisely 2 eigenvalues in the strip < 3^2; < 
h/C , namely and fi{h), where fi{h) is real and of the form 

^i{h) = h{ai{h)e-^^'/'' + a_i(/i)e-2^-^/''), (5-6) 

where aj are real, 

aj{h) ~ ajfi + haj^i + /i — )■ 0, aj^ > 0, 
S, = V{xo)-V{x,). 

As for the problem of return to equilibrium, we obtained the following 
result with F. Herau and M. Hitrik in |40] : 

Theorem 5.3 We make the same assumptions as in Theorem \5.2\ and let Uj 
be the spectral projection associated with the eigenvalue fij, j = 0,1, where 
fiQ = 0, fii = fi{h). Then we have 

= C(l) -.L^^L^ h^ 0. (5.7) 

We have furthermore, uniformly as t > and h ^ 0, 

^-tp/h = iio + e-*^i/'^Hi + C(l)e-*/^, zn C{L\ L^), (5.8) 

where C > is a constant. 
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Actually, as we shall see in the outline of the proofs, these results (as well 
as fl5.5p ) hold for more general classes of supersymmetric operators. 

Very recently we observed with F. Herau and M. Hitrik that we can 
actually treat the case of any (finite) number of local minima. The basic 
observation here is that there is a Hermitian product which becomes a scalar 
product on the space spanned by the Nq lowest eigenvalues, where A^o denotes 
the number of local minima, and for which our operator is formally self- 
adjoint. This makes it possible to apply very much the same methods as for 
the standard Witten complex. 

5.2.3 A partial generalization of | I42] 

Consider on R" {2n is now replaced by n): 



P 




P^+tPi+Po, 



where bj^k,Cj,po are real and smooth. The associated symbols are: 



p(x,0 




Assume, 



P2 >0,Po> 0, 




Assume that 



{x; po{x) = ci(x) = .. = c„(x) = 0} 
is finite = {xi, ...,xn} and put C = {pi, p„,}, pj = {xj, 0). Put 



1 /■^"/^ ^ 

{p)to = 7^ P° exp{tHpJdt, To > fixed. 

^0 J -To/2 
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Here in general we let Ha = a'^ ■ — a'^ ■ denote the Hamilton field of the 
C^-function a = a(x, 

Dynamical assumptions: Near each pj we have {p)to ~ |p ~ PiP ^^^^ 
any compact set disjoint from C we have {p)to ^ (Near infinity this last 

assumption has to be modified slightly and we refer to [39] for the details.) 
The following result from [32] is very close to the main result of [12] and 
generalizes Theorem 15. It 

Theorem 5.4 Under the above assumptions, the spectrum of P is discrete 
in any hand < < Ch and the eigenvalues have asymptotic expansions 
as in ( 15.5)) . 



The linearization of the Hamilton field Hq at pj (for any fixed j) has 
eigenvalues iofc, k = l,..,n with real part ^ 0. Let A_|_ = A+j be the 
unstable manifold through pj for the Hg-fiow. Then A_|_ is Lagrangian and 
of the form ^ = (p'_^_{x) near Xj {(p^ = (p+j), where 



The next result is from |39j : 



Theorem 5.5 Let \j^k{h) he a simple eigenvalue as in / (5.5]) and assume 
there is no other eigenvalue in a disc D{\j^k, h/C) for some C > 0. Then, in 
the sense, the corresponding eigenfunction is of the form e~'^+^^^''^ (a(x; h) + 
0{h°°)) near xj, where a{x; h) is smooth in x with an asymptotic expansion 
in powers of h. Away from a small neighhorhood of Xj it is exponentially 
decreasing. 

The proof of the first theorem uses microfocal weak exponential estimates, 
while the one of the last theorem also uses local exponential estimates. 

5.2.4 Averaging and exponential weights. 

The basic idea of the proof of Theorem 15. 41 is taken from [12], but we reworked 
it in order to allow for non-hypoelliptic operators. We will introduce a weight 
on T*R" of the form 



Put 



?(a;,0 



■p(x, i^) = p2{x, + Pi{x, - Poi^)- 



(j)+{Xj) = 0, 0+(Xj) = 0, 0+(Xj) > 0. 





(5.9) 
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for < e ^ 1. Here J{t) is the odd function given by 

V1-o^<,<., ("0) 

and we choose Pe{p) to be equal to p{p) when dist (p, C) < e, and flatten out 
to ep away from a fixed neighborhood of C in such a way that p^ = 0{e). 
Then 

Hp,llJe= {Pe)To-Pe- (5.11) 

We let e = Ah where A ^ 1 is independent of h. Then the weight 
exp{ip^/h) is uniformly bounded when h 0. Indeed, ipt = 0{h). 

Using Fourier integral operators with complex phase, we can define a 
Hilbert space of functions that are "microlocally 0{e'x^{ipt/h)) in the 
sense". The norm is uniformly equivalent to the one of L^, but the natural 
leading symbol of P, acting in the new space, becomes 

p(exp(2i7^J(p)), pGT^R'^ (5.12) 

which by Taylor expansion has real part ~ Pi{p) + Po(p) + (Pe) ~ Pe- 

Very roughly, the real part of the new symbol is > e away from C and 
behaves like dist(p, C)^ in a A/e-neighborhood of C. This can be used to 
show that the spectrum of P (viewed as an operator on the weighted space) 
in a band < '^z < e/C comes from an -y/e- neighborhood of C. In such a 
neighborhood, we can treat P as an elliptic operator and the spectrum is 
to leading order determined by the quadratic approximation of the dilated 
symbol fl5.12p . This gives Theorem 15.41 

We next turn to the proof of Theorem 15.51 and we work near a point 
Pj = {xj,C,j) G C. Recall that A+ : ^ = (p'^{x) is the unstable manifold for 
the ifg-fiow, where q{x,^) = —p{x,i^). We have q{x, (j)'_^_{x)) = 0. 

In general, if G C°° is real, then P^ := e"^^^ o P o e~^l^ has the symbol 

Pv(^'0 =V2{x,i)-q{x,ii'{x)) +i(g[(x,?/;'(x)) (5.13) 

• As long as g(x, '?/''(a;)) < 0, we have 3f?p^ > and we may hope to 
establish good apriori estimates for P^. 

• This is the case for ip = and for ip = (p^. Using the convexity of g(x, ■), 
we get suitable weights ip with q{x, ip\x)) < 0, equal to (p+^x) near Xj, 
strictly positive away from Xj and constant outside a neighborhood of 
that point. 

• It follows that the eigenfunction in Theorem l5.5l is (roughly) C)(e~'^+(^'^/'*) 
near Xj in the sense. 
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• On the other hand, we have quasi-modes of the form a{x; h)e ^+^^^/'^ 
as in 

• Applying the exponentially weighted estimates, indicated above, to the 
difference of the eigenfunction and the quasi-mode, we then get Theo- 
rem 15.51 

5.2.5 Supersymmetry and the proof of Theorem 15.21 

We review the supersymmetry from [7], [95], see also G. Lebeau [60]. Let 
A{x) : T*R"' — )■ T^R" be linear, invertible and smooth in x. Then we have 
the nondegenerate bilinear form 

{u\v)Aix) = {a''A{x)u\v), u,ve A^^T^R'^, 

and we also write {u\v)a{x) = {u\v)a{x)- 

If u,v are smooth /c-forms with compact support, put 

{u\v)a = j {u{x)\v{x))A{x)dx. 

The formal "adjoint" Q^'* of an operator Q is then given by 

{Qu\v)a = {u\Q^'*v)a. 

Let (j) : R" — )■ R be a smooth Morse function with d"'(f) bounded for 
\a\ > 2 and with |V0| > 1/C for |x| > C. Introduce the Witten-De Rham 
complex: 

= e~h ohdoeh = ^(/i<9x, + dx^(p) o dXj, 
j 

where d denotes exterior differentiation and dxj" left exterior multiplication 
with dxj. The corresponding Laplacian is then: — = d^'*d^ + d(j,d^'*. Its 
restriction to g-forms will be denoted by — Notice that: 

-A5)(e-^/^) = 0. 

Write A = B + C with 5* = B, = -C . -A^ is a second order 
differential operator with scalar principal symbol in the semi-classical sense 
(fa!- ^0) of the form: 

j,k j,k 
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Example. Replace n by 2n, x by (x, y), let 



2 V 7 



Then 

.(0) 



A^^^ = h{ct>'-d,~ct>',-dy) 



y 



When = y^/2 + we recover the KFP operator (15.11) 

The results of Subsubsection 15.2.31 apply, if we make the additional dy- 
namical assumptions there; — A^'* has an asymptotic eigenvalue = o{K) as- 
sociated to the critical point Xj precisely when the index of Xj is equal to 
q (as for the Witten complex and analogous complexes in several complex 
variables). In order to cover the cases g > we also assume that 

A = Const. (5.14) 

The Double well case. Keep the assumption (I5.14p . Assume that is a 
Morse function with {Vcpl > 1/C for |x| > C such that — A^ satisfies the 
extra dynamical conditions of Subsubsection l5T2.3l and having precisely three 
critical points, two local minima U±i and a saddle point Uq of index 1. 

Then — A^"* has precisely 2 eigenvalues: 0, /i that are o{h) while — A^'' has 
precisely one such eigenvalue: fi. (Here we use as in the study of the Witten 
complex, that dtf, and d^'* intertwine our Laplacians in degeree and 1. The 
detailed justification is more complicated however.) e''^^^ is the eigenfunction 
of A^-* corresponding to the eigenvalue 0. Let Sj = 4>{Uq) — 4>{Uj), j = 
±1, and let Dj be the connected component of {x G R"; 0(x) < 4>{Uo)} 
containing Uj in its interior. 

Let E^"^ be the corresponding spectral subspaces so that dimE^^^ = 2, 
dim£'(^) = 1. Truncated versions of the function e~'^^^^^^ can be used as 
approximate eigenfunctions, and we can show: 

Proposition 5.6 E^^^ has a basis ei,e_i, where 

ej = Xj{x)e-TM^'^-<t>(Uj)) ^ ^f^^ E^sense. 

Here, we let Xj G C^{Dj) be equal to 1 on {x E Dj] 0(x) < 4>{Uo) — e}. 
The theorems 15.51 can be adapted to — A^-* and lead to: 
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Proposition 5.7 E^^'^ = Ccq, where 

~ (a; — f/o)^) €0 > is small enough, Qq is an elliptic symbol, Xo ^ 
Co^(R"), Xo = l nearUo. 

Let the matrices of : E^^^ and ci^'* : E^^) with respect 

to the bases {e_i,ei} and {eo} be 

( A_i Ai ) and ^ ^ respectively. 

Using the preceding two results in the spirit of tunneling estimates and com- 
putations of Helffer-Sjostrand ([361 [37]) we can show: 

Proposition 5.8 Put Sj = 0(f/o) — J = ±1- Then we have 

where i±i, are real elliptic symbols of order such that iji* > 0, j = ±1. 

From this we get Theorem 15. 3[ since /i = Al]^A_i + A^Ai. □ 

Thanks to the fact that we have only two local minima, certain simpli- 
fications were possible in the proof. In particular it was sufficent to control 
the exponential decay of general eigenfunctions in some small neighborhood 
of the critical points. 



5.2.6 Return to equilibrium, ideas of the proof of Theorem 15.31 

Keeping the same assumptions, let Ho, Hi be the rank 1 spectral projections 
corresponding to the eigenvalues fiQ := 0, fii := fi of — A^'' and put 11 = 
IIo + Hi. Then e_i^,ei is a basis for 7?.(n) and the restriction of P to this 
range, has the matrix 




(5.15) 
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with the eigenvalues and n = X*_^X_i + A^Ai. A corresponding basis of 
eigenvectors is given by 



vo = ^^(Aie_i - A_ie_i) (5.16) 



vi = ^L(A!iie_i + A*e_i). 



The corresponding dual basis of eigenfunctions of P* is given by 

V* = -L(A*e*_i - Aliel,) (5.17) 
vl = ^(A_ie*_i + Aie*), 

where e*_i,el G TZ{Il*) is the basis that is dual to e_i,ei. It follows that 
Vj,v* = 0(1) in L^, when h 0. 

From this discussion we conclude that Uj = {■\vj)vj, are uniformly bounded 
when h ^ 0. A non-trivial fact, based on the analysis described in Subsub- 
sections I5.2.3[ I5.2.4[ is that after replacing the standard norm and scalar 
product on by certain uniformly equivalent ones, we have 

^{Pu\u)> ^\\uf, VMe7^(l-^), (5.18) 

where 11 is the spectral projection corresponding to the spectrum of P in 
D{0,Bh) for some 5 > 1. 

This can be applied to the study of u(t) := e~*^/'^u(0), where the initial 
state n(0) G is arbitrary: Write 

u{0) = Uou{0) + niu(O) + (1 - n)u(O) =: u° + + u^. (5.19) 

Then 

lk°IUI^^i, Ik^ll < C(1)||m(0)|| (5.20) 

||e-*^/^M^|| < Ce-*/^||M(0)|| (5.21) 

e-'P/\ = e-*'^^/^-, J = 0,1. (5.22) 

Here flCTD follows if we write = {1 - Il)u + (H - Il)u, apply f l538|l to 
the evolution of the first term, and use that the last term is the (bounded) 
spectral projection of n to a finite dimensional spectral subspace of P, for 
which the corresponding eigenvalues all have real part > h/C. □ 
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5.3 Spectral asymptotics in 2 dimensions 

5.3.1 Introduction 

This subsection is mainly based on recent joint works with S. Vii Ngoc and 
M. Hitrik [15] |16], but we shall start by recalling some earlier results that we 
obtained with A. Melin [7T] where we discovered that in the two dimensional 
case one often can have Bohr-Sommerfeld conditions to determine all the 
individual eigenvalues in some region of the spectral plane, provided that 
we have analyticity. In the self-adjoint case such results are known (to the 
author) only in 1 dimension and in very special cases for higher dimensions. 

Subsequently, with M. Hitrik we have studied small perturbations of self- 
adjoint operators. First we studied the case when the classical flow of the 
unperturbed operator is periodic, then also with S. Vii Ngoc we looked at the 
more general case when it is completely integrable, or just when the energy 
surface contains some invariant Diophantine Lagrangian tori. 

5.3.2 Bohr-Sommerfeld rules in two dimensions 

For (pseudo-)differential operators in dimension 1, we often have a Bohr- 
Sommerfeld rule to determine the asymptotic behaviour of the eigenvalues. 
Consider for instance the semi-classical Schrodinger operator 



where we assume that V G C°°(R; R) and V{x) — > +oo, — > oo. Let Eq G 
R be a non-critical value of V such that (for simplicity) {x G R; V{x) < Eq} 
is an interval. Then in some small fixed neighborhood of Eq and for /i > 
small enough, the eigenvalues of P are of the form E = Ek, k E Z, where 



In the non-self-adjoint case we get the same results, provided that '^sV is 
small and V is analytic. The eigenvalues will then be on a curve close to the 
real axis. 

For self-adjoint operators in dimension > 2 it is generally admitted that 
Bohr-Sommerfeld rules do not give all eigenvalues in any fixed domain except 
in certain (completely integrable) cases. Using the KAM theorem one can 
sometimes describe some fraction of the eigenvalues. 

With A. Melin [7T]: we considered an /i-pseudodifferential operator with 
leading symbol p{x,C,) that is bounded and holomorphic in a tubular neigh- 



P 




+ V{x), with symbol p{x,^) = ^'^ + V{x), 




e ■ dx, 9{E; h) ~ ^o(^) + 9i{E)h + ... 
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borhood of R"^ in = Cl x C?. Assume that 



np~^(0) 7^ is connected. 
On R^ we have > 1/C, for \{x,^)\ > C, 



(5.1) 



(5.2) 



for some C > 0, 

d^p{x,^),d'^p{x,^) are hnearly independent for all (x,^) G p~^{0) fl R^. 



(Here the boundedness assumption near oo and (15. 2p can be replaced by a 
suitable ellipticity assumption.) It follows that p~^(0) fl R^ is a compact 
(2-dimensional) surface. 
Also assume that 



"Sufficiently small" here refers to some positive bound that can be defined 
whenever the the other conditions are satisfied uniformly. 

When the Poisson bracket vanishes on p~^{0), this set becomes a La- 
grangian torus, and more generally it is a torus. The following is a complex 
version of the KAM theorem without small divisors (cf T.W. Cherry [T3] . 
J. Moser [73]), 

Theorem 5.9 ( fTT^ ) There exists a smooth 2-dimensional torus T C p~^(0)n 
close to p~^{Q) n R^ such that = and Ij(T) G R, j = 1,2. Here 
Ij(T) := C,-dx are the actions along the two fundamental cycles 7i, 72 C T, 

and a = dC^j A dxj is the complex symplectic (2,0)-form. 

Replacing p hj p — z for 2; in a neighborhood of G C, we get tori 
T{z) depending smoothly on z and a corresponding smooth action function 
I{z) = (/i(r(z)), /2(r(z))), which are important in the Bohr-Sommerfeld rule 
for the eigen-values near in the semi-classical limit /i — )■ 0: 

Theorem 5.10 (fTDjl ) Under the above assumptions, there exists 6q G (|Z)^ 
and e{z; h) ^60 + ei{z)h + e2{z)h^ + .. m C°°(neigh (0, C)), such that for z 
in an h-independent neighborhood of and for h > sufficiently small, we 
have that z is an eigenvalue of P = p{x, hD^) iff 



(5.3) 



is sufficiently small on p ^(0) fl R^. 



(5.4) 



I{z) 



k - e{z- h), for some k G Z\ (BS) 



2'Kh 
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Recently, a similar result was obtained by S. Graffi, C. Villegas Bas |28] . 

An application of this result is that we get all resonances (scattering 
poles) in a fixed neighborhood of G C for —h'^A + V{x) if V is an analytic 
real potential on with a nondegenerate saddle point at x = 0, satisfying 
V{0) = and having {(x, ^) = (0,0)} as its classically trapped set in the 
energy surface {p{x,^) = 0}. 

5.3.3 Diophantine case 

In this and the next subsubsection we describe a result from |3H] and the 
main result of about individual eigenvalues for small perturbations of a 
self-adjoint operator with a completely integrable leading symbol. We start 
with the case when only Diophantine tori play a role. 

Let P^{x, hD; h) on have the leading symbol Pe{x,^) = p{x,^) + 
ieq{x,^) where p, q are real and extend to bounded holomorphic functions on 
a tubular neighborhood of R^. Assume that p fulfills the ellipticity condition 
(15. 2p near infinity and that 

P,=o = Fix, hD) (5.5) 

is self-adjoint. (The conditions near infinity can be modified and we can also 
replace R^ by a compact 2-dimensional analytic manifold.) 

Also, assume that Pe(x, h) depends smoothly on < e < eg with values 
in the space of bounded holomorphic functions in a tubular neighborhood of 
R^, and ~ + ^Pi,e + h^P2,e + when /i — > 0. 

Assume 

p~^(0) is connected and dp ^ Q on. that set. (5.6) 

Assume complete integrability for p\ There exists an analytic real valued 
function / on T*R^ such that Hpf = 0, with the differentials df and dp being 
linearly independent almost everywhere on p~^{0). {Hp = p'^ ■ ^ — p'^ ■ is 
the Hamilton field.) 

Then we have a disjoint union decomposition 

p-^(O) nT*R2 = U ^' ^^•'') 

AeJ 

where A are compact connected sets, invariant under the Hp fiow. We as- 
sume (for simplicity) that J has a natural structure of a graph whose edges 
correspond to families of regular leaves; Lagrangian tori (by the Arnold- 
Mineur-Liouville theorem [22] )• The union of edges J\S possesses a natural 
real analytic structure. 
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Each torus A G J \ S* carries real analytic coordinates Xi,X2 identifying 
A with = R^/27rZ^, so that along A, we have 

d d 

where ai,a2 € R. The rotation number is defined as the ratio w(A) = [ai : 
02] G RP\ and it depends analytically on A G J \ 5. We assume that uj{A) 
is not identically constant on any open edge. 

We say that A G J \ S* is respectively rational, irrational, diophantine if 
ai/a2 has the corresponding property. Diophantine means that there exist 
a > 0, d > such that 

ex 

\{ai,a2)-k\>j^, O^keZ^ (5.9) 

We introduce 

1 /-^/^ 

{q)T = 7f 1° exp{tHp)dt, T > 0, (5.10) 

and consider the compact intervals (5oo(A) C R, A G J, defined by, 

(5oo(A) = [lim inf(g)r, lim sup(g)T]- (5-11) 

A first localization of the spectrum a{P^{x, hD,j.\ h)) ([IS]) is given by 

3(o-(P,) n {z- \^z\ < 5]) C e[inf |J Q^{A) - o(l), sup |J Qoo(A) + o(l)], 

AeJ AeJ 

(5.12) 

when S,e,h 0. 

For each torus A G J \ S', we let {q){A) be the average of with respect 
to the natural smooth measure on A, and assume that the analytic function 
J \ S 3 A (q) (A) is not identically constant on any open edge. 

By combining (15. 8p with the Fourier series representation of q, we see that 
when A is irrational then Qoo{A) = {(g) (A)}, while in the rational case, 

QUA) c (g)(A) + 1], (5-13) 

when uj{A) = — and m G Z, n G N are relatively prime. 

Let Fq G UAej<5oo(A) and assume that there exists a Diophantine torus 
Ad (or finitely many), such that 

(g)(A,) = Fo, rfA(g)(Ad)7^0. (5.14) 

With M. Hitrik and S. Vii Ngoc we obtained: 
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Theorem 5.11 (l4^ ) Assume also that Fq does not belong to Qoo(A) for 
any other A G J. Let < 6 < K < oo. Then 3C > such that for 
h > small enough, and k^^ < e < , the eigenvalues of in the rectangle 
\'Rz\ < h^/C, - e^Fo\ < eh^/C are given by 

P^^\h{k e- h) + 0{h^), k e Z^ 

4 zvr 

Here P^°°\^, e; h) is smooth, real-valued for e = and when h we have 

oo 

P(°°)(e,e;/i) ~ ^(e,6), =p^(0 + ze(g)(0 + O(6^), (5.15) 

corresponding to action angle coordinates. 

In [IS] we also considered applications to small non-self-adjoint pertur- 
bations of the Laplacian on a surface of revolution. Thanks to fl5.13p the 
total measure of the union of all Qoo(A) over the rational tori is finite and 
sometimes small, and we could then show that there are plenty of values Fo, 
fulfilling the assumptions in the theorem. 

With M. Hitrik we are currently studying the distribution of eigenvalues in 
sub-bands that are delimited by two different values "Fo" as in the theorem. 

5.3.4 The case with rational tori 

Let Fo be as in (15.141) but now also allow for the possibility that there is a 
rational torus (or finitely many) A^, such that 

Fo e Qoo{Ar), Fo ^ (g)(A,), (5.16) 

dj,{{q)){Ar) ^ 0, d^{u){Ar)^0. (5.17) 

Assume also that 

Fo ^ goo(A), for all A G J \ {A^, AJ. (5.18) 
With M. Hitrik we showed the following result: 

Theorem 5.12 (fj^) Let 6 > be small and assume that h <^ e < h^~^^ , or 
that the subprincipal symbol of P vanishes and that <^ e < /i3+<^ . Then 
the spectrum of P^ in the rectangle 



[-i,i] + ^e[^o- 
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is the union of two sets: U Er, where the elements of E^ form a distorted 
lattice, given by the Bohr-Sommerfeld rule Ii5.15\) . with horizontal spacing 
X h and vertical spacing x eh. The number of elements i^{Er) of Er is 

NB that X e^+V^^- 

This result can be apphed to the damped wave equation on surfaces of 
revolution. 

5.3.5 Outline of the proofs of Theorem 15.111 and 15.121 

The principal symbol of is p^ = p + ieq + 0{e^). Put 

1 /-^/^ 

{q)t = 7f 1° exp{tHp)dt. 

^ J-T/2 

As in Subsection 15.21 we will use an averaging of the imaginary part of the 
symbol. Let J(t) be the piecewise affine function with support in [— |, |], 
solving 

J'(t) = 5(t)-l[„.,i](t), 

and introduce the weight 

GT{t) = J Ji-^)qoexp{tHp)dt. 

Then HpGr = q— {(i)t, implying 

p, o expiieHc^) =p + ie{q)T + Orie^). (5.19) 

The left hand side of fl5.19p is the principal symbol of the isospectral 
operator e~t'^T(a;,/iDj;)Qp^Qg|-GT(a;,/iD^) g^j^^ under the assumptions of Theorem 

I5.11l resp. I5.12l its imaginary part will not take the value ieF^ on p~^(0) away 
from Arf resp. A^UA^. This means that we have localized the spectral problem 
to a neighborhood of A^ resp. A^ U A^. 

Near A^ we choose action-angle coordinates so that A^ becomes the zero 
section in the cotangent space of the 2-torus, and 

= P(0 + m{x, + 0{e^). (5.20) 

We follow the quantized Birkhoff normal form procedure in the spirit of 
V.F. Lazutkin and Y. Colin de Verdiere [581 [18]: solve first 

i7,G = g(x,0-(g(-,0), (5.21) 
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where the bracket indicates that we take the average over the torus with 
respect to x. Composing with the corresponding complex canonical trans- 
formation, we get the new conjugated symbol 

v{^)+^<q{■.^)) + 0{e^ + ^)■ 

Here the Diophantine condition is of course important. 
Iterating the procedure we get for every A^, 

' v ' 

independent of x 

This procedure can be continued on the operator level, and up to a small 
error we see that is microlocally equivalent to an operator P^{hD^,e; h). 
At least formally. Theorem 15.111 then follows by considering Fourier series 
expansions, but in order to get a full proof we also have take into account that 
we have constructed complex canonical transformations that are quantized 
by Fourier integral operators with complex phase and study the action of 
these operators on suitable exponentially weighted spaces. 

Near we can still use action-angle coordinates as in f l5.20p but the 
homological equation fl5.2ip is no longer solvable. Instead, we use secular 
perturbation theory (cf the book which amounts to making a partial 

Birkhoff reduction. 

After a linear change of x-variables, we may assume that p(^) = ^2 + C^(^^) 
and in order to fix the ideas = ^2 + ^i- Then we can make the averaging 
procedure only in the X2-direction and reduce Pe in ( I5.20p to 

V ' 

independent of X2 , 

where {q)2{xi,^) denotes the average with respect to X2. 

Carrying out the reduction on the operator level, we obtain up to small 
errors an operator Pe{xi, hD^^, hD^^; h) and after passing to Fourier series in 
X2, a family of non-self-adjoint operators on S].^: hD^^, hk; h), k E Z. 

The non-self-adjointness and the corresponding possible wild growth of 
the resolvent makes it hard to go all the way to study individual eigenval- 
ues. However, it can be shown that in the region |^i| ^ e^^^ (inside the 
energy surface p = 0) we can go further and (as near A^) get a sufficiently 
good elimination of the x-dependence. This leads to the conclusion that the 
contributions from a vicinity of A^ to the spectrum of P^ in the rectangle 

\^z\<^, \Qz-eFo\<'-^, 
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come from a neighborhood of A.^. of phase space volume (9(e^/^). 

This explains heuristically why the rational torus will contribute with 
0{e'l'^ j}}}^ eigenvalues in the rectangle. 

The actual proof is more complicated. We use a Grushin problem reduc- 
tion in order to reduce the study near A,, to that of a square matrix of size 
O^e^/"^ /K^). However, even if we avoid the eigenvalues of such a matrix, the 
inverse can only be bounded by 

exp 0(6=^/7^') • (5.22) 

What saves us is that away from A^UA^, we can conjugate the operator with 
exponential weights and show that the resolvent has an "off-diagonal decay" 
like exp(— 1/(C/;.)). This implies that we can confine the growth in f l5.22p to 
a small neighborhood of A^, if 

l_ 

Ch^ h^' 

leading to the assumption e <^ /i^/^ in Theorem 15.121 

Part II 

Non-self-adjoint operators with 
random perturbations. 

6 Zeros of holomorphic functions of exponen- 
tial growth 

We will need a result on the number of zeros in a domain of holomorphic 
functions u{z) = Uh{z) that satisfy an exponential upper bound near the 
boundary of the domain as well as corresponding lower bounds at finitely 
many points, distributed along the boundary. Such a result (related to clas- 
sical results for the zeros of entire functions, cf [61], Chapter III, Section 3, 
Theorem 3) was obtained by Hager [311 [32] under a rather strong regularity 
assumption on the exponent. With Hager [53] we obtained a more general 
result with a logarithmic loss however. Recently I revisited the proofs and 
was able to get a result that includes the earlier ones and allows to eliminate 
such losses, see [91] . 
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Let r d C be an open set and let 7 = 9r be the boundary of F. Let 
r : 7 — j-JO, oo[ be a Lipschitz function of Lipschitz modulus < 1/2: 

\r{x) - r{y)\ <^\x - y\, x,y ej. (6.1) 

We further assume that 7 is Lipschitz in the following precise sense, where r 
enters: 

For every x G 7 there exist new affine coordinates y = (^1,^2) of the 
form y = U{y — x), ?/ G C ~ being the old coordinates, where U = Ux 
is orthogonal, such that the intersection of F and the rectangle '■= {y & 
C; l^il < r(x), 1^2! < Cor{x)} takes the form 

{y e Rx] m > fxim), \yi\ < r{x), } (6.2) 

where fxivi) is uniformly Lipschitz on [— r(x), r(x)], and Co is a fixed con- 
stant, which is larger than the Lipschitz moduli of the functions f^- 

Notice that our assumption (16.21) remains valid if we decrease r. It will 
be convenient to extend the function to all of C, by putting 

r{x) = mf{r{y) -\ — \x — y\). (6.3) 

y^'y 2 

The extended function is also Lipschitz with modulus < |: 

\r{x) - r{y)\ < - y\, x,y eC. 

Notice that ^ 

'"(a;) > -dist (x, 7), (6.4) 

and that 

\y-x\<r{x)^ ^ <r{y)<—^. (6.5) 

Theorem 6.1 Let F (s C 6e simply connnected, and have Lipschitz boundary 
7 with an associated Lipschitz weight r as in l\6.1\) . h6. Put = 
Ux-67-D(x, ar{x)) for any constant a > 0. Let Zj G 7, j G Z/NZ be distributed 
along the boundary in the positively oriented sense such that 

r{z^)/4<\zl,-z'^\<r{z^)/2. 

(Here "4" can be replaced by any fixed constant > 2.) Then if Ci > is 
large enough, depending only on the constant Cq in l\6.^) and if Ci > Cq, 
Zj G D{z'j,r{Zj)/{2Ci)), we have the following: 
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Let (f) be a continuous subharmonic function on ^r/Ci with a distribution 
extension to F U 7r/Ci that will be denoted by the same symbol. Then there 
exists a constant C2 > such that if u = Uh{z) , < h < 1, is a holomorphic 
function on F U 7r/Ci satisfying 

h\n\u\ < (f>{z) on 7r/Ci, (6.6) 

h\n\u{z,)\><P{z,)-e,, for J = 1,2,..., N, (6.7) 

where ej > 0, then the number of zeros of u inT satisfies 

|#(!<-'(o)nr)-^M(r)| < (6^8) 



In — ^\^(dw) 
r{zj) 



Here fi := A0 G P'(F U 7r/Ci) is a positive measure on 7r/Ci so that yu(F) 
and /i(7r/Ci) ore well-defined. Moreover, the constant C2 only depends on Cq 
in 116. 2\} and on Ci. 



We next discuss the elimination of the logarithmic integrals in (16. 8p . 
Using fl6.5l) . we get 

f 2^1 1 / 7 \ L(dz) 

Di.o:-^)JDi.J^) r{z) L{D{zl'-^)) 

\w-z\ L{dz) 
.0, I In — -—\ix{dw) 



where L denotes the Lebesgue measure. Here we use Fubini's theorem and 
the fact that 

/ .(.0, |ln^^|L(dz)<0(l)L(D(.°,^)) 
Jd{zo;-^) r{z) 2Ci 

to conclude that the mean-value of 



\w — z\ 



is 0{l)fi{D{zj , -^)). Thus we can find zj G D{zj, -^^) such that 
^ r I - ~ I 

This leads to the following variant of Theorem 16. 1[ where (16. 8p is simplified 
but where the choice of zj is no more arbitrary. 
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Theorem 6.2 Let F, 7 = dT, r, Zj, Cq, Ci he as in Theorem \ 6.1[ Then 

3 G D{zj,^^) such that if (p, u are as in Theorem \6 . 1\ satisfying h6. 6)) . 
and 

h\n\u{z,)\>(t){z,)-e,, J = 1,2,..., AT, (6.9) 
instead of l[6. 7\), then 



#(«-i(o) n r) - -l-/i(r)| < ^(/i(7.M) + E (6-io) 



2TTh 

Of course, if we already know that 

then we can keep Zj = zj in (16.81) and get fl6.10p . This is the case, if we 
assume that is equivalent to the Lebesgue measure in the following sense: 

, X , on Diz„ #), uniformly for j = 1, 2, iV. 

M/^(.„^)) HDiz,/-^)) 2C,^' 

(6.12) 

Then we get. 

Theorem 6.3 Make the assumptions of Theorem lg.il as well as Ii6.11\) or 
the stronger assumption Ii6.12\) . Then from Ii6.6\) . (6.1), we conclude i\6.10\) . 

In particular, we recover the counting proposition of M. Hager [311 132], where 
was independent of h and of class in a fixed neighborhood of 7. Then 
fi ^ L and if we choose r <^ 1 constant and assume (16. 6p . (16.71) . we get from 

1 r ^ 

mu-\o) n r) - ^/.(r)i <j{r + Yl ^.)- (6-i3) 

1 

Hager had ej = e independent of j, r = y/e, N x e^^/^, so the remainder in 

(ra is o(^). 

We next outline the proof of Theorem 16.11 Using a locally finite covering 
with discs D{x,r{x)) and a subordinated partition of unity, it is standard to 
find a smooth function r(x) satisfying 



^r(x) < r{x) < r{x), \Vr{x)\ < -, d"r{x) = C(f^-I°l). (6.14) 
Cj 2 

From now on, we replace r{x) by r(x) and the drop the tilde. The general 
estimates on r remain valid and we have 

r{x) > — dist (x, 7). 
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Consider the signed distance to 7: 

, . { dist (x,7), X eV . . 

^^^) = |-dist(x,7), xGC\r ^6-^^) 

In the set U^^g^-Rx? we consider the regularized function 

^ — y 

( I \\2 ^i^^) 9{y)L{dy), (6.16) 

where < x ^ C^iD{0,l)), J x{x)L{dx) = 1. Here e > is small and 
we notice that r(x) x r{y), g{y) = 0{r{y)), when xii^ ~ v) / i.^)) 7^ 0- It 
follows that ge{x) = 0{r{x)) and more precisely, since g is Lipschitz, that 

g,{x)-g{x) = 0{er{x)). (6.17) 

Moreover one can show that if (Vg)e denotes the regularization of Vg, ob- 
tained as in (16.161) . then 

V.g,ix)-iVgUx) = Oil) sup (6.18) 

d^g,{x) = 0„((er(x))i-H), |a| > 1. (6.19) 
Let C > be large enough but independent of e. Put 

7c.r = {x e Uy<..iRy- \ge{x)\ < Cer{x)]. (6.20) 

If C > is sufficiently large, then in the coordinates associated to (16.21) . 
7cer takes the form 

fx{yi) <y2< fxivi), \yi\ < r{x), (6.21) 

where are smooth on [— r(x),r(x)] and satisfy 

= 0,{{er{x)Y-'), k > 1, (6.22) 

0</+-/., /x-/-xer(x). (6.23) 
Later, we will fix e > small enough and write 7^ = ^cer and more generally. 

We shall next establish an exponentially weighted estimate for the Dirich- 
let Laplacian in 7^ by adating the general approach of Agmon estimates to 
thin tubes (cf [36], [52]): 
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Proposition 6.4 Let C > be sufficiently large and e > sufficiently small. 
Then if (p e C^(%) and 

\€\ < 7^, (6.24) 



Cr' 



we have 



\\e^Du\\ + < C||re^Au||, u e (H^ H H^)M, (6.25) 

where we use the natural norms. 

Outline of the Proof. Let (p E C'^i'Jr^ R-) put 

-A^ = o (-A) o e-^ = - (^'J^ + .(^'^ oD^ + D^o 

where we make the usual observation that the last term is formally anti-self- 
adjoint. Then for every u G {H"^ fl iJQ)(7r): 

{-A^u\u) = \\D,uf - {{(p',fu\u). (6.26) 

We need an apriori estimate for D^. Let v : % — R" be sufficiently 
smooth. We sometimes consider w as a vector field. Then for u G {H"^ fl 
^o')(7.): 

{Du\vu) — {vu\Du) = i{div {v)u\u). 

Assume div (v) > 0. Recall that if w = Vw, then div (w) = Aw, so it 
suffices to take w strictly subharmonic. Then 

J dw{v)\u\'^dx < 2\\vu\\\\Du\\ < \\Duf + \\vuf, 

which we write 

(div(t;) - \v\'^)\u\'^dx < \\Duf. 



Using this in fl6.26p . we get 

IwDur + 1 (^(div (v) - \v\') - {<p'^nu\'dx < 

||i(-A,)n||||fcn||<i||i(-A,)nf + l||A:nf, 
where k is any positive continuous function on We write this as 

l\\Dur+ y"(i(div(to-ii;r-fc^)-(C)V^<^ii^(-A^)^ir. (6.27) 
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The remaining work is then to see that we can choose v so that 

div(t;) > r"^ \v\ < 0{r^^). (6.28) 

and it turns out that this is possible with 

V = V(e^5/"), (6-29) 

where A > is sufficiently large and g = g^. Then replace v by a small 
multiple and finally choose /c to be a small multiple of 1/r. □ 

If (£ C has smooth boundary, let Gq, Pq denote the Green and the 
Poisson kernels of Q, so that the Dirichlet problem, 

Au = V, u\g^ = f, u,ve C^m, f G C^idQ), 

has the unique solution 



u{x) = / GQ{x,y)v{y)L{dy) + / PQ{x,y)f{y)\dy\. 
Jn Jan 

Recall that -Go > 0, Po > 0. We have 



-Gn{x,y)<C -^\Yi\x-yl (6.30) 

where C > only depends on the diameter of ^2. 
We also have the scahng property: 

G^n(|, |) = GUx, y), x,yem,t> 0. (6.31) 

Moreover, —Gq is an increasing function of Q in the natural sense. Using 
these facts with Proposition 16.41 one can show the following result: 

Proposition 6.5 For all x,y & '~fr (and e > small enough), we have 

- G^^{x,y) < G - -^\n ^^ when \x - y\ < (6.32) 

Zir r[y) C 

1 r^''^'''^ 1 r(y) 
- G^^{x,y) <Gexp{-— —\dt\),when\x-y\>—r, (6.33) 

where it is understood that the integral is evaluated along 7 from vr^(y) & '~f to 
vr^(x) G 7, where 7T^{y) , 7r^(x) denote points in^ with |x— 7r^(x)| = dist (x, 7), 
\y — 7r^(y)| = dist (^,7), and we choose these two points (when they are not 
uniquely defined) and the intermediate segment in such a way that the integral 
is as small as possible. 
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We will also need a lower bound on on suitable subsets of 7,.. For 
e > fixed and sufficiently small, we say that M d 7^ is an elementary piece 
of 7r if 



M C 



• 3?/ G M such that M = y + r{y)M, where M belongs to a bounded set 
of relatively compact subsets of C with smooth boundary. 

In the following, it will be tacitly understood that we choose our elementjiry 
pieces with some uniform control (C fixed and uniform control on the M). 
Using Harnack's inequality one can show: 

Proposition 6.6 If M is an elementary piece in 7^, then 

-G^,(a;,2/)xl + |ln^f|l|, x,y e M. (6.34) 

r{y) 

Let be a continuous subharmonic function defined in some neighborhood 
of 7^. Let 

fi = fi^ = A(f) (6.35) 

be the corresponding locally finite positive measure. 

Let M be a holomorphic function defined in a neighborhood of F U We 
assume that 

h\n\u{z)\ < (f){z), 2; G TV- (6.36) 
Lemma 6.7 Let zq G M, where M is an elementary piece, such that 

h\n \u{zo)\ > 0(20) - e, < e < 1. (6.37) 
Then the number of zeros of u in M is 

<jie + -G^^izo,w)fiidw)). (6.38) 

Proof. Writing as a uniform limit of an increasing sequence of smooth 
functions, we may assume that G C°°. Let 



nu{dz) = ^^27r5(2 — Zj), 
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where Zj are the zeros of u counted with their muhiphcity. We may assume 
that no Zj are situated on d'^r- Then, since A In |n| = Uui 

h\n.\u{z)\ = I G^^{z,w)hnu{dw) + / P^^{z,w)h\n\u{w)\\dw\ (6.39) 
< / G^^{z,w)hnu{dw) + / P^^{z,w)(f){w)\dw\ 

J fr J dfr 

Gy^{z,w)hnu{dw) + (f){z) — / G^^{z,w)fi{dw). 
Putting z = ^0 in (16.391) and using (I6.37p . we get 

G^^{zQ,w)hnu{dw) < e + / -G^^{zQ,w)fL{dw). 

Now ^ 

-G^^{zo,w) >—, w e M, 



and we get (16. 38 p . □ 

Notice that this argument is basically the same as when using Jensen's 
formula to estimate the number of zeros of a holomorphic function in a disc. 
Now we sharpen the assumption (I6.37P and assume 

hln\u{zj)\>(f){zj)-ej, (6.40) 

where Zi,...,Zn- G 7(i_j_-)r are points such that with the cyclic convention 
N+l = l: 

\z,^^-z,\^r{z,), ^^^1- (6.41) 
We also assume that zi, Z2, zn are arranged in such a way that 

Z/NZ I—)- vr-y(zj) runs through the oriented boundary in the positive sense. 

(6.42) 

Let Mj C 7r. be elementary pieces such that 

Zj e Mj, dist {zj, Mk) > when k ^ j, 7^ C UjMj, F = (1 - -^)r. 

(6.43) 

We will also assume for a while that (p is smooth. 
According to Lemma 16. 7^ we have 

#(m-1(0) n Mj) < ^{ej + [ -G^^{zj,w)fx{dw)). (6.44) 
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Consider the harmonic functions on 7fr, 

^(z) = h{\n\u{z)\ + [ -G^^{z,w)nu{dw)), (6.45) 



$(^)=0(2)+/ -G^^{z,w)fx{dw). (6.46) 



Then $(2;) > (f){z) with equahty on (97f. Similarly, \l/(2;) > /iln|M(2;)| with 
equality on d'jr- 

Consider the harmonic function 

H{z) = ^{z) - ^{z), z G -ir. (6.47) 
Then on ^7^, we have by (16.361) that 

E{z) = (piz) - h\n\u{z)\ > 0, 
so by the maximum principle, 

H{z) > 0, on 7^. (6.48) 

By (16.401) . we have 

Hiz,) = ^z,)-^iz,) (6.49) 
= (f){zj) — hln\u{zj)\ 

+ / -G^^{zj,w)iJ,{dw) - / -G^^{zj,w)hnu{dw) 

J 7p J7fr 



< ej+ -G^^{zj,w)iJ,{dw). 
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Harnack's inequality implies that 

H{z) < 0(l)(e, + J -G^^{zj,w)ii{dw)) on M,- n 7^, f = (1 - ^)f. (6.50) 

Now assume that u extends to a holomorphic function in a neighborhood 
of r U 7^. We then would like to evaluate the number of zeros of u in F. 
Using f l6.44p . we first have 

#(n-i(0)n7f)<^$^(e, + ^ -G„Xz,,w)fi{dw)y (6.51) 
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Let X £ C^{T U 7f; [0, 1]) be equal to 1 on F. Of course x will have to 
depend on r but we may assume that for all /c G N, 

V'^x = C(r-^). (6.52) 

We are interested in 

X{z)hnu{dz) = / hln\u{z)\Ax{z)L{dz). (6.53) 



Here we have on 7^ 

h\n\u{z)\ = ^(z) - / -G^-{z,w)hnu{dw) (6.54) 
= ^{z) - H{z) - / -G^^{z,w)hnu{dw) 

= 0(z) + / -G^^{z,w)fi{dw) - H{z) - / -G^^{z,w)hnu{dw) 
= (l){z)+R{z), 
where the last equality defines R{z). 



Inserting this in fl6.53p . we get 

X{z)hnu{dz) = / x{z)lJ'{dz) + / R{z)Ax{z)L{dz). (6.55) 



(Here we also used some extension of to F with /i = A0.) The task is now 
to estimate R{z) and the corresponding integral in fl6.55p . Put 

fij = fi{Mjn-fr). (6.56) 

Using the exponential decay property fl6.33p (equally valid for G^^) we get 
for z e MjD -fr, dist (z, dMu) > r{zj)/0{l), k^j: 

[ -G^^{z,w)fx{dw)< [ -G^^{z,w)ix{dw) + 0{l)J2f^ke'^^'-''[ 

(6.57) 

where \j — k\ denotes the natural distance from j to k in 7a/N7a. Similarly 
from fl6.50p . we get 



H{z) < 0(l)(e, + [ -G^^{zj,w)fi{dw) + ^6-^1^-'^*:), (6.58) 
for z e MjH 7f . 
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This gives the following estimate on the contribution from the first two 
terms in R{z) to the last integral in (^1 



1(1 -G^^{z,w)fi{dw) - H{z)\ Ax{z)L{dz) (6.59) 
= 0(l)^(e,+ / -G,,{z,,w)f,idw)) + Y,' 
+0{1)Y, [ [ -G,^{z,w)f^{dw)\Ax{z)\L{dz). 



/ -G,^\Ax{z)\L{dz) = 0{l), (6.60) 

so fl6.59p leads to 

= 0{1) I ^(7f) + J2^, + Y. I -G^A^v w)fi{dw) 

The contribution from the last term in R{z) (in fl6.54p ) to the last integral 
in (K^ is 

/ / G^^{z,w)hnuidw)Axiz)L{dz). (6.62) 

Here, by using an estimate similar to fl6.57p with fi{dw) replaced by L{dz) 
together with (16.601) . we get 

[ G,^iz,w)iAx)iz)Lidz) = Oil), 

so the expression (16.621) is by (16.511) 

O(/i)#(n-i(0)n7f) (6.63) 

N 

= 0(1) ^(e,+ / i~G,^{z,,w))f^{dw)) 

,•1 J 



N 



0(l)(/i(7,) + ^(e, + / -G^Xzj,w)^l{dw))). 
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This is quite similar to fl6.61p . Using Proposition 16.51 we have 

f \z ■ — w\ 

G^.^{z„w)n{dw) < C(l)( / |ln ' ^ ^ \fi{dw) + i^{Mjn-fr) 

and similarly for the last integral in fl6.63p Using all this in f l6.55p . we get 

X{z)hnu{dz) = J xiz)ij{dz) (6.64) 

+0(l)(/i(7r) + E(^^ + / lln(^^)IMrf^)). 

j J\w-Zj\<r(zj)/C ryZj) 

We replace the smoothness assumption on </> by the assumption that (j) is 
continuous near F and keep (16.401) . Then by regularization, we still get 

Here, we observe that 



|#(«-^(o) nr) - ^1 X{^)hn^{dz)\ < i^{u-\o) 



n7f) 



which can be estimated by means of (16.631) . and combining this with (16.641) . 
we get 

|#(«-i(o)nr)--i-Mr)|< (6.65) 

ZTTtl 



7 The one-dimensional semi-classical case 

In this section we consider a simple model operator in dimension 1 and show 
how random perturbations give rise to Weyl asymptotics in the interior of 
the range of p. We follow rather closely the work of Hager [32] with some 
inputs also from Bordeaux Montrieux [9] and Hager-Sj [33]. Some of the 
general ideas appear perhaps more clearly in this special situation. 

Let P = hD-, + g{x), g e C~(5'^) with symbol p{x,^) = ^ + g{x), and 
assume that '^g has precisely two critical points; a unique maximum and a 
unique minimum. 

Let r2 (E {z G C; min53(7 < '^z < max'^g} be open. Put 

Ps = Ps,. = hD, + g{x) + 6Q^, (7.1) 
Quju{x) = ^ aj^k{uj){u\e'')e\x), 
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where Ci > is sufficiently large, e'^(x) = (27r)^^/^e*'^^, k E 7a, and Uj^k ~ 
A/'(0, 1) are independent complex Gaussian random variables, centered with 
variance 1. is compact, so Ps has discrete spectrum. Let F (s f2 have 
smooth boundary. 

Theorem 7.1 Let n > 5/2 and let > be sufficiently small. Let 5 = S{h) 
satisfy e~'"/^ < 5 < /i*^ and put e = e{h) = h\n{l/S). Then for h > 
small enough, we have with probability > 1 — that the number of 

eigenvalues of Ps in T satisfies 

maiP,) n r) - J-vol(p-^(r))| < Const. ^. (7.2) 
Znri II 

If instead, we let F vary in a set of subsets that satisfy the assumptions 
uniformly, then with probability > l — 0{^) we have fl7.2p uniformly for all 
r in that subset. The remainder of the section is devoted to the (outline of) 
the proof of this result. 

7.1 Preparations for the unperturbed operator 

For z ^ Q, let x+{z),x^{z) G be the solutions of the equation '^g{x) = 
Qz, with ±'^g'{x±) < 0, define ^±{z) by ^± + ^g{x±) = ^z. Then, with 
P± = {x±,C,±), we have 

P(P±) = z, ±^{p,p}{p±) > 0. 
We introduce quasimodes of the form 

ewkb(x) = h-'/^aih)x{x - x+(z))e^*+(^), 

where a{h) ~ ao + hai + .., ao ^ 0, 0+(a;) = Ji^^^^^^iz - g{y))dy, x e 
(neigh (0, R)) and x(x) = 1 in a neighborhood of 0. We can choose a 
depending smoothly on z such that all derivatives with respect to z,^ are 
bounded when /i and ||ewkb|| = 1 where we take the norm over 
]x^{z),x+{z) + 2tx[. We can assume that e„kb is normalized in and 

(P-z)ewkb = O(e-c^). 
Define ^-dependent elliptic self-adjoint operators 

Q = {P- zYiP -z),Q = iP- z){P - zY : L\S^) ^ L\S^), 

with domain P((5), T^{Q) = H'^{S^). They have discrete spectrum C [0, +oo[. 
Using that P — z is Fredholm of index zero, we see that d\mN'{Q) = 
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dim^/{Q). If /i 7^ is an eigenvalue of Q, with the corresponding eigen- 
function e G C°°, then f := {P — z)e is an eigenfunction for Q with the same 
eigenvalue /i. Pursuing this observation, we see that 



Proposition 7.2 T/iere exzsfo a constant C > stic/i i/iai t^ = 0{e 
tf — t^ > h/C for h > small enough. 

Proof. We have Qcwkb = "^'j = 0(e^^/'"'') and since Q is self adjoint 
we deduce that is exponentially small. If Cq denotes the corresponding 
normalized eigenfunction, we see that [P — z)eo =: v with exponentially 
small. Considering this ODE on ]x^{z) — 27r,x_(z)[, we get 



eo 



(x) = C/i-^a(/i)e^^+(^) + Fv{x), 



Fv{x) = ^ I e^('^+(")-<^+(^))t;(y)rf?/, 

where = f^^{z — g{y))dy. We observe that — (p+{y)) > 

on the domain of integration. With some more work Hager showed that 
ll-P'IkcLS.La) = Hence for our particular v, we see that Fv is ex- 

ponentially decaying in L^. Recalling the form of e„kb(3;) we conclude that 
II Co — Cwkbll is exponentially small. 

To show that tf — t^ > h/C, it suffices to show that {Qu\u) > ^||?i||^ 
when M ± Co or in other words, that 

\\u\\ < ^WiP - z)u\\. (7.3) 

If f := (P — z)u, we again have 

u = Ch~-*a{h)eT^^+^''^ + Fv 

for some constant C and the orthogonality requirement on u implies that 

= {l + 0{h^))C+{Fv\eo), 

where {Fv\eo) = 0{h~^)\\v\\, so C = 0{h~^^'^)\\v\\ and we get the desired 
estimate on ||m||. □ 



X 
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7.2 Grushin (Shur, Feschbach, bifurcation) approach 

Let /o be the normalized eigenfunction such that Q/o = ^o/o- observed 
prior to Proposition 17.2^ we get 

(P - z)eo = ao/o, (P - z)* fo = PqCq, ao/^o = 

and combining this with ((P — 2;)eo|/o) = (eo|(P — 2;)*/o), we see that ao = Pq- 
Define R+ : L^{S^) ^ C, P_ = C ^ L^{S^) by 

P+n = (u|eo), R-U- = u„/o. 

Then 



is bijective with the bounded inverse 

£(z) = 



E E+ 

P_ P_4 

Here E = 0{h~^^'^) in — )■ is basically the inverse oi P — z from (/o)"*" 
to (eo)^, E^v = {v\fo), E+V+ = v+eo, P_+ = 0(e-i/(^^)). It is a general 
feature of such auxiliary (Grushin) operators that 

z G a(P) ^ P-+(^) = 0. 

7.3 d-bar equation for 

Proposition 7.3 We have 

d,E.4z) + f{z)E.+ {z)=0, (7.4) 

where 

f{z) = U{z) + f.{z), U{z) = (6tP+)P+, f-{z) = E.d^R.. (7.5) 
Thus, 

a,(e^(^)p_+(^)) = ifd,F{z) = f{z). (7.6) 

Moreover, 

^AF{z) = mdj = li \ - \ ) + 0{1). (7.7) 
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Proof. f l7.4p . f l7.5p follow from the general formula for the differentiation of 
the inverse of an operator, here: 



d^£ + £{ckV)£ = 0. 

Let 11(2;) be the spectral projection of Q: — j- Ccq. It is easy to see that the 
various z and 'z derivatives of Cwkb and Il{z) have at most temperate growth 
in and since eo is the normalization of Il{z)e„]^\, we get the same fact for 
Co and hence for eo — Cwkb- This quantity is also exponentially small in 
and by elementary interpolation estimates for the successive derivatives in 
z,z we get the same conclusion for the higher ^-gradients of Cq — e„kb- 
It follows that 

f+{z) = {eo{z)\d,eo{z)) = {e^khiz)\d,e^i,b{z)) + 0{e~^), 

and the various z, ^-derivatives of the remainder are also exponentially de- 
caying. 

Using that Cwkb behaves like a Gaussian, peaked at the point x+{(), we 
can apply a variant of the method of stationary phase to get 



(ewkb|9.e„kb) = -■^(9,</.+)(x+(z), z) + (7.8) 

where the remainder remains bounded after taking 2;, 2 derivatives. 

Using that (j)+{x+{z), z) = 0, {(f)^y^{x+{z), z) = ^+{z), we get after applying 

dz to the first of these relations, that 

{d,(j)+){x+{z),z) = -^+{z)d,x+{z). 

On the other hand, if we apply dz and to the equation, p{x+{z),^+{z)) = z 
and use that x+ and ^+(2) are real valued we can show that 

{p,p} {p,p} 

Plugging this into f l7.8p . applying and taking real parts, we get the second 
(non-trivial) identity in fl7.7p for the contribution from /_|_. The one from /_ 
can be treated similarly. □ 

Using the expressions for the ^-derivatives of and the analogous 

ones for we have the following easy result relating (17.71) to the sym- 

plectic volume: 
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Proposition 7.4 Writing z = x + iy, we have: 

2 

l{p,P}{p 



dC,+{z) A dx-^{z) = —r — =T7 — 'dy A dx, 



dS,-{z) A dx_{z) = ——f — — TT — -dy A dx, 



2 

:{p,p}{p 

so by 



^AF{z)dy Adx= ^{d^+ A dx+ - A + C(l). (7.9) 
7.4 Adding the random perturbation 

Let X ~ A/c(0,cr^) be a complex Gaussian random variable, meaning that 
X has the probability distribution 

X,{P{duj)) = -^e~^d{^X)d{^X). (7.10) 

Here a > 0. For t < 1/cr^, we have the expectation value 

E(e*l^l^) = (7.11) 
1 — aH 

Bordeaux Montrieux [5] observed that we have the following possibly classical 
result (improving a similar statement in 



Proposition 7.5 There exists Co > such that the following holds: Let 
Xj ~ A/c(0,(tJ), 1 < j < < oo be independent complex Gaussian random 

.2 



variables. Put Si = max cr^. Then for every x > 0, we have 



On \ ^ o X 



p(i:iA-,i=>.)<exp(i±i:4-^), 

Proof. For t < l/(2si), we have 

JV 

P(J]|X,f > x) < E(e*(^l^^l'-^)) = e-*^ J]E(e*l^^l') 

1 

^ 1 

= exp(^^ In 2 t^) ^ 6xp t{Co cr|t — x). 

It then suffices to take t = (2si)~^ □ 
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Recall that 



Q,n(x) = «.-,fc(^)(^|e')e^(x), e\x) = (7.12) 

Since the Hilbert-Schmidt norm of Qoj is given by ||Qw||hs — X] 
we get from the preceding proposition: 

Proposition 7.6 If C > is large enough, then 

C 1 

IIQojIIhs < wi-ih probability > 1 — e~ch^. (7-13) 

Now, we work under the assumption that ||(5w||hs ^ C/h and recall that 
\\Quj\\ < IIQljIIhs- Assume that 

6 < (7.14) 

so that ^ h^^^. Then, by simple perturbation theory we see that 



R+ 

is bijective with the bounded inverse 



= E + Oir) = Oih-^'^) in CiL\ L^) (7.15) 
h 

El = E^ + 0{^) = 0{l) in/:(C,L2) 
Et = E^ + 0{j^) = 0{1) in £(L^ C) 

Ei, = E_+ - 6E_QE+ + 0( 



+ 



As before the eigenvalues of Ps are the zeros of E^^ and we have the d-bar 
equation 

d^Ei^ + f\z)Et^ = 0, 
f{z) = d^R+Ei + Eid^R^ = f{z) + 0{^ ^ 



We can solve d^F^ = (making e^^ Et_^_ holomorphic) with 

^' = ^ + ^(/^) = ^ + ^(]^)l- (7-16) 
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Proposition 7.7 Assume that < t < 1, 5 < /i^/^, 



(5t > e c^o'' 



(5 



(7.17) 



t > 



/i5/2' 



where Cq ^ 1 is fixed. Then with probability > 1 — e cf^ , we have 



(7.18) 



For etJert/ z ^ ^l, we have with probability > 1 — O(t^) — e cP", i/iai 



We only give the main idea of the proof which is to notice that E^Q^Ej^ 
can be written as a sum of independent Gaussian random variables and is 
therefore itself a Gaussian random variable. Applying the standard formula 
for the variance of such a sum we get for the variance: 



where eo(j), /o(j) are the Fourier coefficients of eo, /o- Now we can show 
that the Fourier coefficients are 0{{h/\j\)^) for every > 0, when h\j\ is 
sufficiently large, so if we take Ci (in the definition of Q^) large enough, we 
conclude that = 1 + C(/i°°). 

The remainder of the proof then consists in showing that \E^Q^^Ej^\ is > t 
with probability > 1 — 0{t^) and observing that when this happens, then 
the second term in the expression for E^_j^ in (17.151) is dominant. □ 

Proposition 7.8 Let n > 5/2 and fix eo G]0,1[ sufficiently small. Let 
5 = 5{h) satisfy e"*^"/^ <^ 5 <^ h'^ , and put e = e[h) = hln^. Then with 
probability > 1 — e"^/'-*"'^^^ we have \E^^\ < 1 for all z ^ fl. 
For any z E ^l, we have \Et^\ > e~'"^/^ with probability > 1 — 0{S^/h^). 

This follows from Proposition 17.71 by choosing t such that 



Et4z)\ > 



t6_ 
C 



(7.19) 




(7.20) 



max(-e 



5 




which is possible to do since 



6 
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Under the same assumptions, we also have 

Thus for the holomorphic function u{z) = e^^^^^ Et_^_{z) we have 

• With probabihty > 1 - e-^/^^^'^ we have \u{z)\ < exp{^F{z) + Ce/h) 
for all 2 G fi. 

• For every 2 G fi, we have \u{z)\ > exp(^F{z) — Ce/h) with probability 
> 1 - 0(5V/i5). 

Theorem 17.11 on the Weyl asymptotics of small random perturbations of 
the operator P = hD + g[x) is now a consequence of the following result of 
M. Hager, that we apply with = h^F 

Proposition 7.9 Let F d C have smooth boundary and let (f) he a real valued 
C"^ -function defined in a fixed neighborhood ofT. Let z 1— )■ K) be a family 
of holomorphic functions defined in a fixed neighborhood of T, and let < 
e = e{h) ^ 1. Assume 



\u 



[z; h)\ < exp{j^{(f){z) + e)) for all z in a fixed neighborhood of dV . 



• There exist zi, Z]\f depending on h, with N = N{h) x e such that 
dT C Ui D{zk, y/e) such that \u{zk; h)\ > exp(^(0(2;fc) — e)), 1 < A; < 
N{h). 

Then, the number of zeros of u in T satisfies 

|#(«"i(0) nF) [ A<j>{z)dxdy\ < C^. 

This is essentially a special case of Theorem 16. 3^ but we outline the simple 
and direct proof of Hager in the next subsection. 



7.5 Proof of Proposition 17. 9L an outline 



Define (pj{z) by i(pj{z) = (f>{zj) + 2dz4>{zj){z — Zj). Then 

(t){z) = 3f?(20,(2)) + Rj{z), Rj{z) = 0{{z - ZjY) 

<f>'^{z) = '^dz^{z) + 0{{z-zj)). 
Consider the holomorphic function 
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Then \vj{z; h)\ < eiW^)-3**</'i W) = e^^^ < e^, when z - Zj = 0{^/e), while 

\Vj[Zj; h)\ > e . 

In a A/e-neighborhood of 2j we put v = Vj and make the change of variables 
w = {z — Zj) / y/e, v{w) = v{z), so that 

\v{w)\ < e^'l^ on £'(0,2), \v{^)\ > e~^'/\ 



Using Jensen's formula we see that the number of zeros wi, wat of w in 
-D(0,3/2) (repeated with their multiplicity) is 0{e/h). Factorize: 



TV 

v[w) = e^^"'^ ~ ^k) 



1 



Using the maximum principle and a suitably chosen disc of radius between 
4/3 and 3/2, and then also Harnack's inequality we can follow a standard 
procedure to show that 

^g{w), g'{w) = 0{e/h) in D(0,6/5). 

Using finally that dT is covered by the discs D{zj,y/e) and using the 
above representation of u in each disc, we can show that the number of zeros 
of u{-] h) in r is equal to 

^ f ''^'^''^z = 3fJ-\ / -d,<P{z)dz + oA 



2m Jqy u{z) 2iTh Jqy i h 



27ih 



J A(j){x)dxdy + 0{^). □ 



8 The mult i- dimensional semi-classical case 
8.1 Introduction 

In this section we consider general semi-classical operators with multiplicative 
random perturbations. We follow [HH |89] which make use of the work |33j . 
The use of Theorem 16.21 rather than the corresponding weaker result in [53] 
led us to improved remainder estimates in comparison with [55] . 

Let X be a compact smooth manifold on which we choose a positive 
density of integration so that the scalar product on L'^{X) is well-defined. 
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On X we consider an /i-differential operator P which in local coordinates 
takes the form, 

P= <^a{x;h){hDr, (8.1) 

|Q|<m 

where we use standard multiindex notation and let D = D,,. = l-^. We 

t ax 

assume that the coefficients Cq, are uniformly bounded in C°° for h g]0, ho], 
< /lo ^ 1- (We will also discuss the case when we only have some Sobolev 
space control of ao{x).) Assume 

a„(x;/i) = a°(x) + C(/i) in (8.2) 
aa{x; h) = aa{x) is independent of /i for |a| = m. 

Notice that this assumption is invariant under changes of local coordinates. 

Also assume that P is elliptic in the classical sense, uniformly with respect 
to h: 

for some positive constant C, where 

\a\=m 

is invariantly defined as a function on T*X. It follows that pm{T*X) is a 
closed cone in C and we assume that 

pUT*X) ^ C. (8.5) 

If ^ C \pm(r*^), we see that A^o ^ S(p) if A > 1 is sufficiently large and 
fixed, where S(p) := p{T*X) and p is the semiclassical principal symbol 

p{x,0= Yl «°(^)^" (8.6) 

\a\<m 

Actually, fl8.5p can be replaced by the weaker condition that S(p) 7^ C. 

Standard elliptic theory and analytic Fredholm theory now show that if 
we consider P as an unbounded operator: L'^{X) — )■ L'^{X) with domain 
'D[P) = H^{X) (the Sobolev space of order m), then P has purely discrete 
spectrum and each eigenvalue has finite algebraic multiplicity. 

We will need the symmetry assumption 

P* = rPr, (8.7) 
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where P* denotes the formal complex adjoint of P in L'^{X, dx), and dx is the 
fixed smooth positive density of integration and F is the antilinear operator 
of complex conjugation; Tu = u. Notice that this assumption implies that 

p(x,-0=p(x,0, (8.8) 

and conversely, if p fulfills flS.Sp . then we get (18. 7p if we replace P by \{P + 
rP*r), which has the same semi-classical principal symbol p. Actually, (18. 7p 
can be formulated more simply by saying that P is symmetric for the bilinear 
form u{x)v{x)dx. 

Let vXt) ■■= vol ({p e T*X; \p{p) - < t}). For k e]0, 1], ^ G C, we 
consider the property that 

V,{t) = O(t^), < t < 1. (8.9) 

Since r i— )■ p{x, r^) is a polynomial of degree m in r with non- vanishing 
leading coefficient, we see that (18. 9 p holds with k = l/(2m). 
The random potential will be of the form 

qUx) = ^ ak{u;)ek{x), \a\cD < R, (8.10) 

0<tik<L 

where is the orthonormal basis of eigenf unctions of h'^R, where R is an 
/i-independent positive elliptic 2nd order operator on X with smooth coeffi- 
cients. Moreover, h^Reu = fJ^l^k, fJ'k > and we may assume for simplicity 
that the /i^ form a (non-strictly) increasing sequence. We choose L = L{h), 
R = R{h) in the interval 

h-^^ <^L<Ch~^\ M> (8.11) 

s ^ e 

^/,-(t+^)M+.-f < ^ < ch-^, M>|^_^+(^ + e)M, 

for some e g]0, s — |[, s > |, so by Weyl's law for the large eigenvalues of 
elliptic self-adjoint operators, the dimension D is of the order of magnitude 
(L/h)^. We introduce the small parameter 6 = ro/i^^"*^", < tq < Vh, where 

' — ' T) 

Ni:=M + sM+-. (8.12) 
2 

The randomly perturbed operator is 

P5 = P + 5/i^ig^=:P + 5g^. (8.13) 
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The random variables aj{uj) will have a joint probability distribution 

P{da) = C(/i)e*("'^)L(rfa), (8.14) 

where for some A^4 > 0, 

|V,$| = 0{h-^^), (8.15) 

L{da) is the Lebesgue measure and we use the standard i"^ norm on C'^. 
{C{h) is the normahzing constant, assuring that the probability of Bcd (0, R) 
is equal to 1.) 

We also need the parameter 

eo{h) = (/i- + /i"ln^)(ln- + (ln^)2) (8.16) 
n tq h 

and assume that tq = tqQi) is not too small, so that e^Qi) is small. Let 
r2 (E C be open, simply connected, not entirely contained in S(p). The main 
result of this section is: 

Theorem 8.1 Under the assumptions above, letT ^VL have smooth bound- 
ary, let K g]0, 1] be the parameter in /i8.10\) . /i8.11\) . /i8.16\) and assume that 
lis. 9\) holds uniformly for z in a neighborhood of dV . Then there exists a con- 
stant C > such that for > r > 0, ?> CeQ{h) we have with probability 

> 1 ^ ^^oC') 8.17 

that: 

maiPs) n r) - ^^vol ip~\m < (8.18) 

^(^ + vol(j9-i(9r + D(0,r)))^. 

Here i^{a{Ps) fl F) denotes the number of eigenvalues of Ps in T, counted 
with their algebraic multiplicity. 

Actually, the theorem holds for the slightly more general operators, ob- 
tained by replacing P by Pq = P + ^oih^Qi + 12)^ where 
||?2||i/'' < 1, < 5o < ^- Here, is the standard Sobolev space and 

is the same space with the natural semiclassical /i-dependent norm. See 
Subsection 18.31 This allows us in principle to consider more general random 
perturbations and will be used in Section [Hi 
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We also have a result valid simultaneously for a family C of domains 
r C r2 satisfying the assumptions of Theorem 18.11 uniformly in the natural 
sense: With a probability 

>1 \ ) o\ I 8.19 

^2J^n+raax{n{M+l),NA+M) ^ 

the estimate f lS.lSp holds simultaneously for all T EC. 

Remark 8.2 If k > 1/2, then volp"^((9r + £'(0,r)) = 0{r'^''-^), where the 
exponent 2ft; — 1 is > 0. More generally, if 

vol {j)-\dV + D(0, r))) = 0(r°), 

for some a e]0, 1], then we can choose r = and obtain that that the right 
hand side in (18.181) is (!?(l)/i~"'?«+T showing that we have Weyl asymptotics. 
Notice here that if z is not a critical value of p, in the sense that d^p{p) and 
Q'p(p) are independent whenever p{p) = z, then (18. 9 p holds with k, = 1. 

In the proof we replace the zero counting proposition from [33] by the 
stronger Theorem 16.21 leading to an improved remainder estimate. It may 
be possible (though we have not yet checked the details) to replace the right 
hand side in (I8.18P by 

^voi(p-i(9r + D(o,/ii-^))), 

for any fixed e > 0, and also to let F be /i-dependent of a suitable Lipschitz 
class as in section [61 

Remark 8.3 When R has real coefficients, we may assume that the eigen- 
functions ej are real. Then (cf Remark 8.3 in [SB]) we may restrict a in (l8.1Up 
to be in R^ so that is real, still with |a| < R, and change C{h) in (I8.14p 
so that P becomes a probability measure on B^d^O, R). Then Theorem 18.11 
remains valid. 

Remark 8.4 The assumption (18. 7p cannot be completely eliminated. In- 
deed, let P = hDx + g{x) on T = R/(27rZ) where g is smooth and complex 
valued. Then (cf Hager [31]) the spectrum of P is contained in the line 
^z = ^g{x)dx / {2ti). This line will vary only very little under small mul- 
tiplicative perturbations of P so Theorem 18.11 cannot hold in this case. On 
the other hand, for other classes of perturbations, like the ones in Section [7] 
or in [33], the symmetry assumption can be dropped. 

In the remainder of this section, we shall outline the proof of Theorem 18.11 
following [891 EH]. 
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8.2 Semiclassical Sobolev spaces and multiplication 

We let if^(R"') C 5'(R"), s G R, denote the semiclassical Sobolev space of 
order s equipped with the norm where the norms are the ones in 

L^, i"^ or the corresponding operator norms if nothing else is indicated. Here 
{hD) = (1 + {hOyy^'^. In [SB] we recalled the following result: 

Proposition 8.5 Let s > n/2. Then there exists a constant C = C{s) such 
that for all u,v e Hf^{R"), we have u G L°°(R"), uv G Hf^{W) and 

\\u\\lo. <Ch-''''/^\\u\\Hf^, (8.20) 

\\uv\\Hf^<Ch-''/^\\u\\HsJv\\Hs^. (8.21) 

We cover X by finitely many coordinate neighborhoods Xi,...,Xp and 
for each Xj, we let xi, ...,Xn denote the corresponding local coordinates on 
Xj. Let < Xj ^ C^{Xj) have the property that J2iXj > on X. Define 
H^{X) to be the space of all u G V{X) such that 

p 

\H\ ■■= Yl WxAhDYXjuf < oo. (8.22) 
1 

It is standard to show that this definition does not depend on the choice of the 
coordinate neighborhoods or on Xj- With different choices of these quantities 
we get norms in (18.221) which are uniformly equivalent when /i — > 0. In fact, 
this follows from the /i-pseudodifferential calculus on manifolds with symbols 
in the Hormander space S'j^g, that we quickly reviewed in the appendix in [SB]. 
See also [BH], Section 4. An equivalent definition of H^^X) is the following: 
Let 

h^R = J2ihD^,yrj,k{x)hD,, (8.23) 

be a non-negative elliptic operator with smooth coefficients on X, where 
the star indicates that we take the adjoint with respect to the fixed posi- 
tive smooth density on X. Then h^R is essentially self-adjoint with domain 
H'^{X), so (1 + h'^Ry/'^ : — )■ is a closed densely defined operator for 
s G R, which is bounded precisely when s < 0. Standard methods allow 
to show that (1 + h'^Ry/'^ is an /i-pseudodifferential operator with symbol 
in SIq and semiclassical principal symbol given by (1 + r(x,^))'^/^, where 

r{x,^) = k'rj^k{x)^j^k is the semiclassical principal symbol of h'^R. See 
the appendix in [88] . The /i-pseudodifferential calculus gives for every s G R: 
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Proposition 8.6 H^{X) is the space of all u G V'{X) such that {l+h'^Ry/'^u G 
and the norm \\u\\h'> is equivalent to ||(1 + h'^RY^'^u\\, uniformly when 

Remark 8.7 From the first definition we see tliat Proposition 18.51 remains 
valid if we replace R*^ by a compact n-dimensional manifold X. 

Of course, H^{X) coincides with the standard Sobolev space H^{X) and 
the norms are equivalent for each fixed value of /i, but not uniformly with 
respect to h. The following variant of Proposition 18.51 will be useful when 
studying the high energy limit in Section O 

Proposition 8.8 Let s > n/2. Then there exists a constant C = Cs > 
such that 

\\uv\\h^^ < C\\u\\hs\\v\\h^^, \fu G H'{YC), v G HliYC). (8.24) 
The result remains valid if we replace R" hy X. 

The proof is straight forward. We work in local coordinates and make a 
Fourier transform. Then we have to estimate convolutions in certain weighted 
spaces. See [89] for the details. 

8.3 i7^-perturbations and eigenf unctions 

Let S'^{T*X) = S'^oiT*X), S'^iU x R") = S'^.^iU x R") denote the classical 
Hormander symbol spaces, where U C R" is open. The condition (18. 5p 
implies that the closure of the image of p is not equal to the whole complex 
plane and (as in [321 EH] we can find p G S'^{T*X) which is equal to p outside 
any given fixed neighborhood of p~^{VL) such that p — z is non- vanishing, for 
any z & VL. Let P = P + Opf^(j)—p), where Opf^(j)—p) denotes any reasonable 
quantization of (j) — p){x, h^). (See for instance the appendix in [88].) Then 
P — z : HJ^{X) — H^{X) has a uniformly bounded inverse for z ^ Q and 
h > small enough. Now (see for instance [331 [SS]) the eigenvalues of P in 
Q, counted with their algebraic multiplicity, coincide with the zeros of the 
function z ^ det((P^- z)-\P - z)) = det(l - (P - z)-^{P - P)). Notice 
here that (P — z)"^{P — P) is of trace class so the determiant is well-defined 

Fix s > n/2 and consider the perturbed operator 

Ps = P + Sih^q^ + q^) = P + SiQ, + Q^) = P + SQ, (8.25) 
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where qj G if'^(X), 

\\qi\\Hi < 1, IMhs < 1, < 5< 1. (8.26) 

According to Propositions 18.51 18. 8[ Q = 0{1) : — H'^ and hence by 
duahty and interpolation, 

Q = C(l) -.H^-^H^, -s<a< s. (8.27) 

Again, the spectrum of Ps in Q is discrete and coincides with the set of 
zeros of 

det((n - z)-\Ps - z)) = det(l - (Ps - zy\P - P)), (8.28) 

where := P + 5Q. Here (P - z)'^ = 0{1) : ^ for a in the same 
range and by an easy perturbation argument, we get the same conclusion for 

{Ps-z)-'. 
Put 

P,,, := {Ps - zY\Ps -z) = l-{Ps- z)-\P - P) =: 1 - Ks,.^ (8.29) 

Ss,, := Pl^Ps,, = 1 - (K^,. + - Kl^Ks,;) =: 1 - L^,,. (8.30) 
Clearly, 

Ks,,,Ls,, = 0{l):HY ^Hl (8.31) 

For < a < 1/2, let tTq, = l[Q^a]{.Ss^z)- Then using some simple functional 
calculus we showed in [SB], that 

vr, = 0{l) : ^ (8-32) 

We have the corresponding result for Ps — z. Let 

Ss = {P6-zr{Ps-z) (8.33) 

be defined as the Friedrichs extension from C°°{X) with quadratic form do- 
main HJ^{X). For < a < C(l), we now put tTq = l[o,a](5'5). Then as in [88] . 
we see that this new spectral projection also fulfils fl8.32p . for < a <^ 1. 

8.4 Some functional and pseudodifFerential calculus 

Let P be of the form (18.11) and let p in (18. 6p be the corresponding semi- 
classical principal symbol. Assume classical ellipticity as in (18.31) and let 
2; G C be fixed throughout this subsection. Let 

S = {P-z)*{P-z), (8.34) 
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viewed as the self-adjoint Friedrichs extension from C°°. Later on we will 
also consider a different choice of S, namely 

S = P*P,, where P, = {P - zy\P - z) (8.35) 

and P is defined prior to (18.251) . The main goal is to make a trace class 
study of x{\S) when 0</i<a^l, C^(R). With the second choice 
of S, we shall also study lndet(S' + ax(^'S')), when x ^ 0, x(0) > 0. The 
main step will be to get enough information about the resolvent (w — ^S*)"^ 
for w = 0{1), '^w 7^ and then apply the Cauchy-Riemann-Green-Stokes 
formula _ 

xi-S) = -- [ ^^iw - -Sr'Lidw), (8.36) 
a IT J ow a 

where x ^ C'(f'(C) is an almost holomorphic extension of Xi so that 

|| = C(|53w;r). (8.37) 
cfw 

Thanks to fl8.37p we can work in symbol classes with some temparate but 
otherwise unspecified growth in l/IQ'wI. 
Let 

s = \p-z\^ (8.38) 

be the semiclassical principal symbol of S in (18.341) . A basic weight function 
in our calculus will be 

l + S 



A:=(^) , (8.39) 



satisfying y/a < A < L 

As a preparation and motivation for the calculus, we first consider symbol 
properties of 1 + ^ and its powers. 

Proposition 8.9 For every choice of local coordinates x on X, let (x,^) 
denote the corresponding canonical coordinates on T*X . Then for all £ G R, 
a,f3 & N", we have uniformly in ^ and locally uniformly in x: 

afaf (1 + -Y = 0(1)(1 + -)^A-l"l-l^l(e)-l^l. (8.40) 

The proof ([SH]) is straight forward and the same can be said about the proof 
of 

Proposition 8.10 (fS^) Let w vary in some bounded subset of C. For all 
£ G R, 5, /3 G N", there exists J G N, such that 

&Aw - -Y = C(l)(l + -)'A^I"l^l''l(e)~'^'|^>w^|"'', (8.41) 

uniformly in ^ and locally uniformly in x. 
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We now define our new symbol spaces. 

Definition 8.11 Let m(x,^) be a weight function of the form m(x,^) = 
(^)'=A^ We say that the family a = e C°°(T*X), w e D{0,C), belongs to 
S'A(m) if for all 5, /3 G N" there exists J G N such that 

d^d^a = O{l)m{x,0A-^^^-^^^{0'^'^^\^w\-^. (8.42) 

Here, as in Proposition I8.10[ it is understood that the estimate is ex- 
pressed in canonical coordinates and is locally uniform in x and uniform in 
Notice that the set of estimates fl8.42p is invariant under changes of local 
coordinates in X. 

Let U C X be a coordinate neighborhood that we shall view as a subset of 
R" in the natural way. Let a G Sa{T*U, m) be a symbol as in Definition 18.111 
so that fl8.42p holds uniformly in ^ and locally uniformly in x. For fixed values 
of a, w the symbol a belongs to Siq(T*U), so the classical /i-quantization 

Au = Op,{a)u{x) = J I e^^^-yy^a{x, r^; h)u{y)dydT^ (8.43) 

is a well-defined operator C^{U) C°°{U), £'{U) V'{U). In order to 
develop our rudimentary calculus on X we need a pseudolocal property for 
the distribution kernel KA{x,y), whose proof is also routine (see [89]). 

Proposition 8.12 For all a, (3 e N", G N, there exists M G N such that 

d^d^KAix, y) = 0{h^\'^wr''), (8.44) 

locally uniformly on U x U \ diag(f/ x U). 

This means that if G C^{U) have disjoint supports, then for every 
iV G N, there exists M G N such that Mijj : H-^{W) H^{W) with 
norm 0{h^\'^w\~'^^), and this leads to a simple way of introducing pseudo- 
differential operators on X: Let Ui, ...,Us be coordinate neighborhoods that 
cover X. Let Xj ^ C^{Uj) form a partition of unity and let Xj ^ C^{Uj) 
satisfy Xj ~< Xj i^i the sense that Xj is equal to 1 near supp(xj)- Let a = 
[Oil , . . . , Oig ), where aj G Sx{fn). Then we quantize a by the formula: 

s 

^ = 5Zx,°0p;,(a,)ox,. (8.45) 

1 

This is not an invariant quantization procedure but it will suffice for our 
purposes. 

Using integration by parts and stationary phase we can study the compo- 
sition to the left with non-exotic pseudodifferential operators and we obtain 
the following result for a coordinate neighborhood: 
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Proposition 8.13 (;89^). Let A = Op;,(a), a E 5i,o(mi), B = Oph{b), b E 
Sa{'^2) o-nd assume that b has uniformly compact support in x. Then AoB = 
Op^(c), where c belongs to S!^{mim2) and has the asymptotic expansion 

in the sense that for every N ElSi, 

7,1/31 

c= J2 -^d^^a{x,OD^Xx,0 + rN{x,^;h), 

where r^ESf^i^^h^). 

We have a parametrix construction for w — j^S, still with S as in fl8.34p . 
Let us first work in a coordinate neighborhood U, viewed as an open set in 
R". Then for every G N we can construct a symbol 

such that on the symbol level 

_ ^S)#E^ = 1 + r^, e 5a((-^)^+^), (8.47) 

E^- is a holomorphic function of w, for |^| > C, (8.48) 

where C is independent of A^. 

Now we return to the manifold situation and denote by E^^\ rj^^ the 
corresponding symbols on T*Uj, constructed above. Denote the operators 
by the same symbols, and put on the operator level: 

s 

En = Y,XjE^nX,. (8.49) 

i=i 

with Xj, Xj ill (18.451) . Then 

{w-^S)E^^, = l-Y^}-[S,x,]Ef^rX,+Y.X,4\, (8.50) 
=: 1 + -Rj^^ + b!-^ 



l + B 



■N- 
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Proposition 18.121 implies that for every A^, there exists an M such that the 



trace class norm of R)^ satisfies 



(1) ^^^^ 



As for the trace class norm of -Rjy , we can combine standard facts about 



\R)^'\\tr<0{h"\Qw\-'''). (8.51 

5(2) 

such norms for pseudodifferential operators and scaling to get 

||i?iv||tr < 1 1' (^-^y dxdC (8.52) 

The contribution to this expression from the region where A > 1/C is 
C(/t^-'^)|3u;|-*^(^). 

The volume growth assumption (18. 9|) . that we now assume for our fixed 
z, says that 

V{t) := vol ({p e T*X; s < t}) = Oit"), < t < 1, (8.53) 
for < K < 1. Using this and (I8.52p one can show that 





N 






la. 





\\RN\\tr < 0{l)h-V ( - ) (8.54) 
From f l8.50p . we get 

{w - -Sy^ = En-1 -{w- -S)-^Rn. 
a a 

Write ^ ^ 

More precisely we do this for each E^^_^ in (18.491) . Then quantize and plug 
this into flOHD: 



X{-S) = -- [ ^Op,{^-^)L{dw)-- [ ^F^.,L{dw) (8.55) 



a 



-- f ^{w - -Sy^RNL(dw) =: I + II + III. 

71 J OW a 

Here by definition, 



w / — ' \w 
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with the coordinate dependent quantization appearing to the right. 
After some further estimates we get 

As at the last estimate in the proof of Proposition 4.4 in [33] we see that this 
quantity is 0{a'^h~^) and more generally, 

For II, one can show, using the fact that the symbol is holomorphic in w 
for large ^, that 

a 

It is also clear that 

a' /-h'" 



Summing up our estimates, we get the following result: 
Proposition 8.14 Let x G C^(R). For 0<h<a<l,we have 

||x(^5)|k. = 0(l)^, (8.57) 

trx(-5) = -^ II xC^)dxd^ + 0(^-). (8.58) 
a (zTTAi)" J J a nP- a 

Remark 8.15 Using simple h-pseudodifferential calculus (for instance as in 
the appendix of fgg| /. we see that if we redefine S as in ^8.35^) . then in each 
local coordinate chart, S = Opfj^{S), where S = s mod Sifi{h{^)~^) and s is 
now redefined as 



2 



The discussion goes through without any changes (now with m = 0) and we 
still have Proposition \8. 14\ with the new choice of S, s. 

In both cases it follows from ([g.57| j that the number N[a) of eigenvalues 
of S in the interval [0,a] satisfies 

0{a^/h''). (8.60) 
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In the remainder of this subsection, we choose S, s as in fl8.35p . fl8.59p . In 
this case we notice that is a trace class perturbation of the identity, whose 
symbol is 1 + 0{h°° / {^)°°) and similarly for all its derivatives, in a region 
\i\ > Const. 

Let < X G e'er ([0 5 with x(0) > and let ao > be small and fixed. 
Using standard pseudodifferential calculus in the spirit of [7D], we get 

\ndetiS + aoxi—S)) = -^i[[ Hs + aoxi—s))dxd^ + Oih)). (8.61) 
ao (27r/i)" J J ao 

Extend x to be an element of C^(R; C) in such a way that t + x(t) ^ 
for all t G R. As in [53], we use that 

JtHE + tx{j)) = \^P{j), (8.62) 

where 

so that ip G C(j"(R). By standard functional calculus for self-adjoint opera- 
tors, we have 

|lndet(5 + tx(f))=tri^(^). (8.64) 
Using fl8.58p . we then get ior t > a > h > 0: 

Integrating this from t = ao to t = a and using fl8.6ip . (18. 62 p . leads to 

Proposition 8.16 If < x G C^([0, oo[), x(0) > 0, we have uniformly for 
< /i < a < 1 

lndet(^ + ax(-5)) = , \, ( [ [ In s{x, ^dxd^ + O {a'' In a)). (8.65) 
a (27r/i)" J J 

Here the remainder term can be replaced by 0{a'^) when k < 1 and by 0{a + 
kin a) when k = 1. 

Notice that (I8.65P imphes the upper bound, 

lndetP;P, <-4— (// \nis)dxd^ + Oia'^\n-)). (8.66) 
(27r/i)" J J a 
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We next consider P^^^ = {P5 - z^^iPs - z) = 1 - Ks,^ with Ps = P + 5Q, 
Ps = P + SQ as in Subsection 18.31 and under the assumptions fl8.28p . (18.301) . 
Put 

Ss,. = Pl.Ps,. = 1 - Ks,. - + K^Ks,., 
where Kg^z is given by (18.301) . so that 

\\KsJ<0{l), \\KsJtr< im - z)-'\\\\P - PI, < Oih-'^). 

Here || • ||tr denotes the trace class norm, and we refer for instance to [23] for 
the standard estimate on the trace class norm of an /i-pseudodifferential op- 
erator with compactly supported symbol, that we used for the last estimate. 
Write ks,, = f,Ks,.. Then 

ks,z = -{z - Psr'Qiz - Psr\p - p), 

so 

\\ksJ<0{\\Q\\), \\ksJtr<0{\\Q\\h-^). 

It follows that 

\\SsA<0{\\Q\\), \\Ss4tr<0{\\Q\\h~^). 

Let = N{a, 6) denote the number of singular values of Ps^z in the 
interval [0, ^/a[ for h a <^ 1. Assume 

6 < 0{h). (8.67) 

Then \\Ss,z - 5o,^|| < 0{h) and from frail) we get 

Ar(a,5) = C(a'^/i-"). (8.68) 

Let la(t) = max(a,t). For < e < 1, let C°°(R+) 3 la,e > be equal 
to t outside a small neighborhood of t = and converge to 1^ uniformly 
when e — 0. For any fixed e > 0, we put f(t) = lQ,e(t) for t > and extend 
/ to R in such a way that /(t) = t + g{t), g e C^(R). Let f{t) = t + g{t) 
be an almost holomorphic extension of / with g G C^(C). Then we have: 

f{Ss,z) = Ss-^j{w- Ss,z)''dg{^)L{dw). 

Differentiating with respect to 6 one can show the identity 

^lndet/(S',,,) = tTif{Ss,,y'f'iSs,z)Ss,.). 
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Now we can choose / = lo,,^ such that \f'{t)\ < 1 for t > 0. Then we get the 
estimate 

lndet(l„,,(55,,)) = tT{l^,,{Ss,.rH'^^,{Ss,.)Ss,.) 
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/-/-)/ II 'S'(5,z ||tr ■ 

a ' 

WQW 



0(1) 



Since Indet 1q(S'5^z) = hm^-s^o In det la,e(5'5^z), we can integrate the above 
estimate, pass to the hmit and obtain 

Indet = Indet l„(^o,.) + ' 



With some more work, we also get 



Indet = , \, ( / / lns(x,Ot^a;rf^+CKlna)+C(^^)). (8.69) 

8.5 Grushin problems 

Let P : H — 7- "H be a bounded operator, where 7/ is a complex separable 
Hilbert space. Following the standard definitions (see [27]) we define the 
singular values of P to be the decreasing sequence si(P) > S2{P) > ... of 
eigenvalues of the self-adjoint operator (P*P)i/2 as long as these eigenvalues 
lie above the supremum of the essential spectrum. If there are only finitely 
many such eigenvalues, Si(P), Sfc(P) then we define Sfc+i(P) = s,fc_|„2(-P) = 
... to be the supremum of the essential spectrum of (P*P)^/^. When dim?/ = 
M < oo our sequence is finite (by definition); si > S2 > ... > sm, otherwise 
it is infinite. Using that if P*Pu = s'jU, then PP*{Pu) = s'jPu and similarly 
with P and P* permuted, we see that Sj{P*) = Sj{P). Strictly speaking, 
P*P : Af{P)^ Af{P)^ and PP* : Af{P*)^ Af{P*)^ are unitarily 
equivalent via the map P[P*P)-^I'^ ; M(P)-^ — )■ M(P*)-^ and its inverse 
P*{PP*)-y^ : Af{P*)^ Af{P)^. (To check this, notice that the relation 
P{P*P) = {PP*)P on 7V(P)^ implies P{P*PY = {PP*YP on the same 
space for every a G R.) 

In the case when P is a Fredholm operator of index 0, it will be convenient 
to introduce the increasing sequence < ti(P) < t2{.P) < ••• consisting first 
of all eigenvalues of (^P*PY/^ below the infimum of the essential spectrum and 
then, if there are only finitely many such eigenvalues, we repeat indefinitely 
that infimum. (The length of the resulting sequence is the dimension oiT-L.) 
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When dim'H = M < oo, we have tj{P) = SM+i-jiP)- Again, we have 
tj{P*) = tj{P) (as reviewed in [33] )• Moreover, in the case when P has a 
bounded inverse, we see that 

SAP-') = J^y (8.70) 

Let P be a Fredholm operator of index 0. Let 1 < < oo and let 
R+ : Ti — 7- C^, -R_ : — 7- H be bounded operators. Assume that 

V=( !^ V X ^ -H X (8.71) 



is bijective with a bounded inverse 

'-{It) <«-^^' 

Recall (for instance from [91]) that P has a bounded inverse precisely 
when E |_ has, and when this happens we have the relations, 

P-^ = E- E+EzlE^, Ezl = -R+P-^R-. (8.73) 

Recall ([27]) that if A,B are bounded operators, then we have the general 
estimates of Ky Fan, 

Sn+k-M + B)< Sr,{A) + Sk{B), (8.74) 

Sn+k~i{AB) < Sn{A)skiB), (8.75) 
in particular for /c = 1, we get 

Sn{AB) < \\A\\s^{B), Sn{AB) < Sn{A)\\Bl SniA + B)< s„(A) + 

Applying this to the second part of fl8.73p . we get 

SkiEzl) < \\R4\\R+\\sk{P-'), l<k<N 

implying 

tkiP) < \\R^ II \\R+\ME_+), l<k<N. (8.76) 

By a perturbation argument, we see that this holds also in the case when P, 

E ^ are non-invertible. 

Similarly from the first part of ( 18.73 p . we get 

Sk{P^')<\\E\\ + \\E4\\E4s,{Ezl), 
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leading to 

'''^> ^ ii^iiM^-:f;ikiiii^-ir 

Again this can be extended to the non-necessarily invertible case by means 
of small perturbations. 

Next, we recall from [33] a natural construction of an associated Grushin 
problem to a given operator. Let Pq : Ti ^ H he a. Fredholm operator of 
index as above. Assume that the first singular values ti(Po) ^ t2{Po) ^ 
... < t7v(-fo) correspond to discrete eigenvalues of PqPq and assume that 
tN+i{Po) is strictly positive. In the following we sometimes write tj instead 
of tj{Po) for short. 

Recall that are the first eigenvalues both for PqPq and PqPq. Let 
ei, and /i, /at be corresponding orthonormal systems of eigenvectors 
of PqPq and PqPq respectively. They can be chosen so that 

Poe,=tjf„ P^f,=tje,. (8.78) 

Define R+ : L'^ ^ and i?_ : ^ by 

N 

RMj) = (w|e,), R-U- = J2^-U)fj- (8-79) 



As in [33], the Grushin problem 



PqU + -R_M_ 

R+u = v^, 



(8.80) 



has a unique solution (u, n_) G x for every {v, v^) G x C^, given 

by 



where 



= J]i;+(j)e„ E%{j) = {v\f,), (8.82) 
1 

i5;°^ = -diag(t,), ll^°ll<i^- 

i?" can be viewed as the inverse of Pq as an operator from the orthogonal 
space (ei,e2, ...,eAr)-^ to (/i, /2, /tv)"^- 

We notice that in this case, the norms of R^ and R_ are equal to 1, 
so fl8T6|l tells us that tfc(Po) < 4(^°+) for 1 < A; < A^, but of course the 
expression for E^_^ in (18.821) implies equality. 
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Let Q G C{H,'H) and put Ps = Pq — SQ (where we sometimes put a 
minus sign in front of the perturbation for notational convenience). We are 
particularly interested in the case when Q = Q^^u = q^iU is the operator of 
multiplication with a random function q^. Here 5 > is a small parameter. 
Choose R± as in fl8.79p . Then \i 5 < t^^i and \\Q\\ < 1, the perturbed 
Grushin problem 

Psu + R^u^ = V, 
R^u = v.. 



is well posed and has the solution 

u = 1 

U-= Et+ Et+v+, 



u = Eh) + 

n<5 I 7-i<5 



is obtained from by 



1 

oo 



1 

oo 



(8.83) 



where 

^'=( fs ) (8.85) 



Using the Neumann series, we get 

Et^ = + 6E°_QEl + 5^E°_QE^QEl + 6^E'^_Q{E^QfEl + ... (8.87) 
We also get 

oo 



El = El + J2 S^E^QfEl (8.89) 



Ei=E'^^ + J2 iQE'^)^- (8.90) 
1 

The leading perturbation in E^^ is SM, where M = E^QE^ : ^ 
has the matrix 

M{uj)j,k = (Qcklfj), (8.91) 
which in the multiplicative case reduces to 



M{u)j,k= I q{x)ek{x)fj{x)dx. (8.92) 
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Put To = tAr+i(Po) and recall the assumption 

ligil < 1. (8.93) 

Then, if 5 < ro/2, the new Grushin problem is well posed with an inverse 
given in flSysHjl -f lOOjl . We get 



\EH < -^^IEH < — , \\Ei\\ < -^-^ < 2, (8.94) 



To 

52 1 52 



11^1+ - (^°+ + SE'^^QEDW < ^— < 2 — . (8.95) 

ro 1 - - To 

Using this in (glE>, (IgTrj) together with the fact that tk{Et+) < 2to, we get 

^^^^ < h{Ps) < h{Ei^). (8.96) 

o 

Remark 8.17 under suitable assumptions, the preceding discussion can be 
extended to the case of unbounded operators. This turns out to be the case 
for our elliptic operator P5. 



We next collect some facts from [33] . The first result follows from Section 
2 in that paper. 

Proposition 8.18 Let P : H — ?■ H be bounded and assume that P — 1 is of 
trace class, so that P is Fredholm of index 0. Let R^, R-,V,S = be as 
in (8/71), (8/72). Then V is also a trace class perturbation of the identity 
operator and 

det P = det P det E^+. (8.97) 

Now consider the operator P^ = Pq^z in (18.291) for 2; G fi, and keep the 
assumption fl8.53p . 
Define 

P<5,z R~,5 
R+,5 

as in ([HIZHD-dHIHDD, so that V = Vo- As in (5.10) in [33] we have 



detVsl^ = a^^ detlo,{S5,z), 2 In | det P^] = Indet 1„(55,^) + iVln -, 



a 

(8.98) 



where la(t) = max(a,t), t > 0, which with fl8.69p and the bound 
0{a'^h~'^) gives 



In I det P^l = jjln \pz\dxd^ + 0{a^ In i + ^IIQII))- (8. 



99) 
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8.6 Singular values and determinants for certain ma- 
trices associated to 5 potentials 

Lemma 8.19 (JM TO Let ei, G (7°(X) and put 

( ei(x) \ 



62 (a:) 
V e7v(a;) / 



X e X. 



Let L C he a linear subspace of dimension M — 1, for some 1 < M < N . 
Then there exists x & X such that 

dist (^(x), Lf > ^;^j^tr ((1 - txl)Sx), (8.100) 

where Ex = {{ej\ek)L^(x))i<j,k<N and ttl is the orthogonal projection from 
onto L. 

Proof. Let i^i, z/at be an orthonormal basis in such that L is spanned 
by ui, vm~i (and equal to when M = 1). Then by direct calculations, 



N 

i dist {'t{x),Lfdx = ^{Sxi^iWi) = tr ((1 - 7Tl)Sxj 



It then suffices to estimate the integral from above by 

vol (X) sup dist {~^{x), LY, 

and we can find an x G X satisfying flS.lOOp . □ 
If we make the assumption that 

61, is an orthonormal family in L^(X), (8.101) 

then Sx = ^ and (18.1001) simplifies to 

msixdist(-^(x),Lf > . ,,,, . 8.102 

x£X VOl(X) 

In the general case, let < ei < 62 < ■■■ < ej\f denote the eigenvalues of 
Sx- Then, using the mini- max principle, one can show that 

inf tT{{l-7rL)£x) =ei + e2 + ...+eN-M+i=: Em. (8.103) 

dim L=A/— 1 
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Now, we can use the lemma to choose successively ai,...,aN G X such 
that 

ll^(ai)f > 



vol (X) ' 



dist(^(a2),C^(ai))2 > 



vol (X) ' 



dist(^(aM),C^(ai)©... ©C^(aA/-i))' > 



vol (X) ' 



Let 1/1,1/2, ■■■,1'N be the Gram-Schmidt orthonormalization of the basis 
(oat), so that 



"e^(aM) = CM^Mraod(jyi, ...,pm~i), where |cm| > (^ y^fl'^^ ^ • (8-104) 

Consider the N x N matrix E = ("e^(ai) ~^{a2) ... ~^{aj\i)) where ~^{aj) 
are viewed as columns. Expressing these vectors in the basis z/i,...,z/jv will 
not change the absolute value of the determinant and E now becomes an 
upper triangular matrix with diagonal entries Ci, ...,CAr. Hence 



and fl8.1U4p implies that 



det^l = |ci ■ ... ■ ctvI, (8.105) 



Let /i, /2, /at be a second family of continuous functions on X. Define 
M = ^ by 



N 



Mu = ^(M|7(a^))-e^(a^), u E C^. (8.107) 
1 

Then 

M = EoF*, (8.108) 

where 

F=(t(ai)...?(a;v)). (8.109) 

Now, we assume 

fj = e„yj. (8.110) 
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Then F* =^E, so 

M = Eo^E. CS.lir 



We get from flSJOHD . fISTTTD . that 



vol (X)^ 

Under the assumption flS.lOip . this simphfies to 



detM|> -^'f;-;f^" . (8.112) 



Using that 

N 



detM\ = Y[sj< st^s^''''^^ < sf , where Sj = Sj{M), (8.114) 



we get 

Proposition 8.20 Under the above assumptions, 



- vol(X) ' 



^'^^■(n(;;;^)j ■ ("IS) 

8.7 Singular values of matrices associated to suitable 
admissible potentials 

We let P, P,p,phe as in the introduction to this section. We also choose e^, 
/ifc, D = D{h), L = L{h) as in and around fl830D . fl8ll]) . 

Definition 8.21 An admissible potential is a potential of the form 

q{x) = J2 «fcefc(x), a E C^. (8.117) 

o<^lk<L 

We shall approximate (5-potentials in H^'^ with admissible ones and then 
apply the results of the preceding subsection. As in the introduction we let 
s > n/2, < e < s - n/2. 
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Proposition 8.22 Let ae X. Then 3a e , r e H^"" such that 

W = Yl "fcefc + r(x), (8.118) 

where 

Mh-^ < C,,L~(^-t-^)/i-t, (8.119) 

h ' 

iJ2 l^/cD^ < (L)^+^( (A^;c)-^(^+^)|«.r)^ < CLt+^/.-t. (8.120) 

The proof uses that 5a G if with norm 0{h~"'^'^) and the fact that 

for any s G R, 

oo 

II ^ttfcefcll^j ^ y^(/^fc)^''|afc|^- 
1 

It then suffices to truncate the expansion of 5a in the basis ei, e2, .... 

Let Ps be as in IK^ and assume ([H2SD, flOU]) . Let 7^(7^«) = Cci © ... © 
Cejv be as in one of the two cases of subsection l8.3[ By the mini-max principle 
and standard spectral asymptotics (see [23]), we know that = 0{h~"') and 
if we want to use the assumption (18.91) we even have = (9((max(a, h))'^h~^) 
bv fl8.68l) . For the moment we shall only use that is bounded by a negative 
power of h. Let a = (ai, a^v) G and put 

N 

ga(a^) = 5^5(x-a,), (8.121) 
1 

Mg^;,-fc = j qa{x)ek{x)ej{x)dx, l<j,k<N. (8.122) 

Now IK32\i implies that \\Y. >^kek\\Hj^ < C(||A||£2) so flCTD and the fact that 
= C(l)Ar;i-"/2, imply that for all A,/i e C^, 

(Mg^A,/i) = j ga(x)(^Afcefc)(^/ijej)c?a; 
= C(l)iV/i^"||A||||^|| 

and hence 

si{Mj = ||MgJU(civ,ciV) = (9(l)Ar/i-". (8.123) 

We now choose a so that fl8.115p . fl8.116p hold. 

The Cj form and orthonormal system, so = 1 and 

Ej = N - J + 1. (8.124) 
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Then, choosing the aj as in Subsection 18.61 (18.1 ISp gives the lower bound 

where the last identity follows from Stirling's formula. 
Rewriting fl8.116p as 



N „ \ JV-fc+i 



and using (18.123^ . we get 



Sk> , ...... . (-) (iV!)^. (8.126) 



CN-k+i (vol {X))N~k 
Summing up, we get 



+1 



Proposition 8.23 We can find ai, cat G X such that if Qa = 

and Mq^.j^k = / (la{x)ek{x)ej{x)dx , then the singular values Si > 82 > ... > 

sn ofMgl, satisfy [HJM) . [SlJm and [SAM) . 

We next approximate with an admissible potential by applying Propo- 
sition to each (5-function in q^. 

qa = q + r, q= ^ a^efc, (8.127) 

where 

WqWh-^ < Ch-^N, (8.128) 

h 

Mh-^ <C,L~'^'~^~'^h-^N, (8.129) 

(Yl l"'^!')^ - CL'^^'h-^N. (8.130) 
Below, we shall have = 0{h'^~^) so if we choose L as in (18. lip , we get 

and q satisfies flgrTU]) . flSlT]) . We get 

||M,|| < C,L^('^t-^)/i-"Ar. (8.131) 
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For the admissible potential q in fl8.127p . we thus obtain from ( I8.126p . 

mm- 



CiefctT(voi(x))^wi 

Similarly, from (]8123|1 . flSTTHT]) we get for L > 1: 



(ATI) iTOT _ CeL^(^"t-^)/i-"Ar. 

(8.132) 



|Mg|| < CNh- 



(8.133) 



(8.134) 



Using Proposition 18. 6[ we get for every e > 0, 

WqWh^^ < 0{l)NL'+^+'h-^, Ve > 0. 
Summing up, we have obtained 

Proposition 8.24 Fix s > n/2 and Ps as in |07| j, [KM) . [KM) and let 



TTa, ei,...,e]\f be as in one of the two cases in Subsection \8.3[ Choose the 
h-dependent parameter L with 1 ^ L < 0{h~^'^) for some fixed Nq > 0. 
Then we can find an admissible potential q as in ^8.127^ (different from the 
one in \8.21\j , \8.2t^] ) such that the matrix Mg, defined by 



M, 



qCkCjdx, 



satisfies ^8.132\) . \8.13'J\) . Moreover the H^-norm of q satisfies \8.134\j. 



Notice also that if we choose R with real coefficients, then we can choose 
q real- valued. 



8.8 Lower bounds on the small singular values for suit- 
able perturbations 

In this subsection, we fix a z G ^2. We will use Proposition 18.241 iteratively to 
construct a special admissible perturbation for which we have nice lower 
bounds on the small singular values of P5 — 2, that will lead to similar bounds 
for the ones of Ps^z and to a lower bound on | det -P<5,2|. 
We will need the symmetry assumption (18. 7p : 

P* = VoPoV, (8.135) 

This property remains unchanged if we add a multiplication operator to P. 
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As in the introduction, we let 

V,{t) = vol ({p G T*X- \p{p) - < t}), (8.136) 
and assume (for our fixed value of z) that (18.91) holds: 

V,{t) = 0{t^), < t < 1, (8.137) 
for some k, e]0, 1]. Proposition 18.141 gives: 

Proposition 8.25 Assume ^8.131^ and recall Remark \8. 11\ For < h <^ 

a ^ 1, the number N{a) of eigenvalues of (P — z)*{P — z) in [0, a] satisfies 

N{a) = Oia^'h-''). (8.138) 
Let e>0,s>| + ebe fixed as in the introduction and consider 

Po = P + Hh'^qi + ^2°), with < 5o < hlUi. hlWns < 1. (8.139) 

From the mini-max principle, we see that Proposition 18.251 still applies after 
replacing P hj Pq. 

Choose To g]0, ((7/1)2] and let = 0{h'^~^) be the number of singular 
values of Pq — ^; < ti(Po — z) < ... < tisi{Po — z) < tq in the interval [0, to[. 
As in the introduction we put 

Ni=M + sM + -, (8.140) 

where M, M are the parameters in (18. lip . Fix 6 g]0, ^[ and recall that A^ is 
determined by the property tj^{PQ ~ z) < tq < t7v+i(Po — z). Fix eo > 0. 

Proposition 8.26 a) If q is an admissible potential as in i\8.10\} . 118.11]} . we 

have ^ 

Iklloo < Ch-'^/^WqllHi < Ch-^\ (8.141) 

b) If N is sufficiently large, there exists such an admissible potential q, such 
that if 

P, = P,+ ^JfLq=:P, + 5Q, 5 = ^/.^^+" 
(so that \\Q\\ < 1 then 

t,{Ps-z)>t,iPo-z)- °^ > (i-__)t,(Po-^), u>N, (8.142) 
tu{Ps -z)> roh^\ [N -eN] + l<u <N. (8.143) 
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Here, we put 

N2 = 2{Ni + n)+ eo, (8.144) 

and we let [a] = max(Zn] — 00, a]) denote the integer part of the real number 
a. When N = 0{1), we have the same result provided that we replace ^^KJJ^ 
by 

tNiPs) > roh^^. (8.145) 



Proof. The part a) follows from Subsection 18. 2[ the definition of admissible 
potentials in the introduction and from the definition of Ni in fl8.140p . (See 
also fl8.134p .) We shall therefore concentrate on the proof of b). 

Let ei,...,eAr be an orthonormal family of eigenfunctions corresponding 
to ty{PQ — -z), SO that 

(Po - zy{P^ - z)ej = it,{Po - z)fej. (8.146) 

Using (18. 6p -v^ fl8.135p . we see that a corresponding family of eigenfunctions 
of (P — z){P — z)* is given by 

f,=Te,. (8.147) 

/i,...,/Ar and fi,.-.,fN are orthonormal families that span the same space 
P/v. Let En be the span of ei, cn- We then know that 

(Po - z) : En ^ Fn and (Pq - z)* : Fn ^ En (8.148) 

have the same singular values 0<ti<t2<---<tN- 
Define R+ : C^, P_ : ^ by 

N 

RMj) = (w|e,), R-U^ = Y,U-U)fr (8-149) 

1 

Then 

Pq — Z P„ \ . ^ / D \ X. r^N > r2 X, r^N 



P = Q J : P(Po) X C^' ^L'x (8.150) 

has a bounded inverse 

P _ ( E P+ 

The singular values of P_+ are given by tj{E^^) = tj{Po — z), I < j < N, 
or equivalently by Sj(P_+) = tN+i-jiPo — z), for 1 < j < N. 

We will apply Subsection 18.51 and recall that N is assumed to be suffi- 
ciently large and that 6 has been fixed in ]0, l/4[. (The case of bounded 
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will be treated later.) Let N2 be given in fl8.144l) . Since z is fixed it will also 
be notationally convenient to assume that z = 0. 

Case 1. Sj{E_+) > Toh^\ for 1 < j < iV - [(1 - e)N]. Then we get the 
desired conclusion with q = 0, Ps = Pq. 

Case 2. 

Sj{E^+) < Toh^^ for some j such that 1 < J < iV - [(1 - 0)N]. (8.151) 

Recall that for the special admissible potential q in (18.1271) . we have 
(Kim . For k < N/2, we have - A; + 1 > N/2, so 

k-1 

< 1, 



N-k + 1 
and fl8132|l gives 

^.(M,)>^(iV!)^-aL-(^-5-)-. 

By Stirling's formula, we have (A^!)^ > A^/Const, so for 1 < k < N/2, we 
obtain with a new constant C > 0: 

Here, we recall from Proposition 18.251 (which also applies to Pq) that = 
0{h'^~"') and choose L so that 

i.e. so that (in agreement with (18. lip ) 

L^h^^^'. (8.152) 

We then get 



Sk{Mg) >—, l<k<—, (8.153) 



h"" , N 

for a new constant C > 0. 

From (I8.133|) and the fact that A^ = C(/i^-") we get 

si(Mg) < IIM^II < CA^/i"" < Ch^-'^''. (8.154) 

In addition to the lower bound fl8.152p we assume as in (18. lip (in all 
cases) that 

L < Ch-^, for some M > ^ • (8-155) 



s- 2 -e 
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As we saw after fl8.130p . q is indeed an admissible potential as in flS.lUp . 
dHUD, so that by flCTTj) 



q\\ oo <Ch- 2 \\q\\Hs^<Ch~^\ (8.156) 



Put 



Ps = Po + -^q = Po + 6Q, Q = — g, ligil < 1. (8.157) 

c c 

Then, if 5 < ro/2, we can replace Pq by Ps in fl8.150p and we still have 
a well-posed problem as in Subsection 18.51 with = Q as above. Here 
E^_QEl = h^'Mq/C so according to fl8J[53D . we have with a new constant 
C 

Sk{5E'_QEl) > 1 ^ ^ ^ y (8.158) 

Playing with the general estimate f l8.74p . we get 

and for a sum of three operators 

s,{A + B + C)> s,+k+i_2{A) - Sk{B) - se{C). 
We apply this to Et_^ in fl8.95p and get 

s.iE'_^) > s,+k-i{5EtQEl) - Sk{E^_^) - 2-. (8.159) 

Here we use fl8J[5T]) with j = k = N - [{I - e)N] as well as f l8J[58|) . to get 
{oiv<N- [(1 - e)N] 

sAEU) > - roh'^^ - 2-. (8.160) 

O To 

Recall that 6 < \. 
Choose 

5 = ^ro/i^^+", (8.161) 

where (the new constant) C > is sufficiently large. 

Then, with a new constant C > 0, we get (for /i > small enough) 

Su{Et+) > 1 < < iV - [(1 - e)N], (8.162) 
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implying 

s^{Ei^) > 8ro/^^^ 1 < u < N - [{1 - e)N]. (8.163) 
For the corresponding operator Pg, we have for u > N: 



UiPs) > U{Po) -S = U{Pq) - 



C 



Since ty{P) > tq in this case, we get fl8.142l) . 
From mm and ^M), we get fimg]) . 

When iV = we still get fl8160|) with z/ = 1 and this leads to (ISlMHl) . 

□ 

The construction can now be iterated. Assume that ^ 1 and re- 
place (Po,iV,ro) by (P5, [(1 - ^)iV], tq/i^^) =: (pW, iV^^), r^^^) and keep on, 
using the same values for the exponents Ni,N2. Then we get a sequence 
{P^'^\ JV^s)^ Tq''^), /c = 0, 1, k{N), where the last value k{N) is determined 
by the fact that N^''^^^'' is of the order of magnitude of a large constant. 
Moreover, 

t,{P^^^) > iV^'^) <u< iV^^-i), (8.164) 
^^(p(fc+i)) > ^^(pW) _ ToJL , > atW, (8.165) 

ATC^ + I) = [(1 -^)iVW], (g_;Lg7) 

p(0) ^ p^ ^(0) ^ ^(0) ^ 

Here, 

p(fc+l) ^ p{fc) + ^(fe+l)g(fc+l) ^ p(fe) ^ 

Notice that A^'^'^'^ decays exponentially fast with k: 

Nik) < (1 _ o)''N, (8.168) 
so we get the condition on k that (1 — 6)'^N > C ^ 1 which gives, 

In^ 

k < - — ^. (8.169) 
1-6* 



We also have 



^0 



(fc) 



ro(/^'^^)^ (8.170) 
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For u > N, we iterate fl8.165p . to get 

uNi+n 

UiP^"^) > t.(P)-ro-^(l + /i^^ + /i2^^ + ...) (8.171) 

uNi+n 

For 1 ^ < A^, let £ = ^(A^) be the unique value for which A^^^-* < < 
iV(^-i), so that 

tu{P^'^)>4'\ (8.172) 

by (18164]) . lik> i, we get 

l.Ni+n 

U{P^'^) > K{P^'^) - ri'^0{—^). (8.173) 

In ^ 

The iteration above works until we reach a value k = = 0{ , 'f ) 

for which A^^*-'") = (9(1). After that, we continue the iteration further by 
decreasing A^*^'^^ by one unit at each step. 

Summing up the discussion so far, we have obtained 

Proposition 8.27 Let (P, z) satisfy the assumptions in the beginning of this 
subsection and choose Pq as in ^8AM)- Let s > , < e < s-^ , M > 

Ni = M+sM+|, N2 = 2{Ni+n)+eo, where eo > 0. Let L be an h-dependent 
parameter satisfying 



n — 3n 

-M 



h^^^ < L < Ch^''\ (8.174) 

Let < To < \/h and let N'^^^ = (9(/i^^") be the number of singular values of 
Pq — z in [0, Tq[. Let Q < 9 < j and let N{9) 1 be sufficiently large. Define 
N^'^\ 1 < k < ki iteratively in the following way. As long as N^''^ > N(9), 
we put A^^'^+i) = [(1 — 9)N^''^]. Let ko > be the last k value we get in this 
way. For k > ko put A^C^+i) = A^^*-') — 1, until we reach the value ki for which 

Put Tq^"* = TqH'^^^, 1 < k < ki + 1. Then there exists an admissible 
potential q = qh{x) as in li8.1(J\) . liS.ll]) . satisfying li8.iyU\) . I^ETI^^, so that, 

Mh^ < 0(l)/i"^^+t, IkllL- < 0(l)/^-^^ 

such that ifPs = Po + ^ro/i2^^+"g = Pq + SQ, 6= ^h^'+"To, Q = h^^q, we 
have the following estimates on the singular values of Ps — z: 

• Ifu> A^(o), we have t^{Ps - z) > {1 - ^^%{Po - z). 
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• //M^-) <u< M^-i), l<k<ki, thenU{Ps-z) > {l-0{h^'+''))T^''\ 

• Finally, for u = N^''^^ = 1, we have hiPs-z) > (1 - C(/i^i+"))r^'''+^\ 

Now it is possible to pass from the Grushin problem for P^ — 2; to a suitable 
one for P^^ and follow up with estimates on the singular values (cf fl8.77p ) 
and obtain: 



Proposition 8.28 Proposition 8.27 remains valid if we replace Ps — z there 
with P^ z- 

Taking a suitable Grushin problem for P^^z and using Proposition 18. 181 we 
can show when Tq = \/h: 

Proposition 8.29 For the special admissible perturbation Ps in the propo- 
sitions \KW^ \8.28\ we have 

ln|detP5,^|> (8.175) 

We also have the upper bound 

|detE_+| < <exp(CiV(°)), 

which leads to 

In I det Ps,z\ < (^jjln \pz\dxd^ + O + Z^" In ^ 

(8.176) 

Notice that this bound is more general, it only depends on the fact that 
the perturbation of P is of the form 5Q with 5 = tqH^^^^/C and with 

\\Q\\=o{i). 

When To < \/h we keep the same Grushin problem as before and notice 
that the singular values of P-+ that are < tq, obey the estimates in Propo- 
sition |8271 Their contribution to ln|detii^ |_| can still be estimated from 

below. The contribution from the singular values of E y. that are > tq can 

be estimated from below by — (!?(/i''^" ln(l/ro)) and this leads to the conclu- 
sion that Proposition \8.2fA remains valid when < tq < \/h. The same holds 
for the upper bound ^8.170^. 
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8.9 Estimating the probability that det Et_^ is small 

In this subsection we keep the assumptions on (P, z) of the beginning of 
Subsection 18.81 and choose Pq as in (18.1391) . We consider general Ps of the 
form 

P5 = Po + SQ, 5Q = 5 = ^/i^^+Vo, (8.177) 

where q is an admissible potential as in fl8.10p . (18. lip . Notice that D : = 
#{^; A^fc < L} satisfies: 



D < 0{L^h~'^) < 0{h^^'), Ns := n{M + 1) 
With R as in (I8.10p . we allow a to vary in the ball 

\a\cD < 2R = 0{h~^). 



(8.178) 



(8.179) 



(Our probability measure will be supported in B(jd{0,R) but we will need 
to work in a larger ball.) 

We consider the holomorphic function 



F{a) = (detP5,,)exp( 



{2Trhy 



In \pz\dxd^). 



Then by (18.1761) . we have 

ln|F(a)| < eo(/i)/i"", |a| < 2R, 



(8.180) 



(8.1811 



and for one particular value a = with \a^\ < ^R, corresponding to the 
special potential in Proposition 18.271 



ln|F(a°)| > -eo{h)h- 



(8.182) 



where eo{h) is given in (I8.16p . 

Let e C'^ with = R and consider the holomorphic function of one 



complex variable 



f{w) =F{a° + wa^). 



(8.183) 



We will mainly consider this function for w in the disc determined by the 
condition la" + wa^\ < R: 



w + 



\R' R 



< 1 



a 



R 



+ 



P"' P' 



=: r, 



0; 



(8.184) 



whose radius is between ^ and 1. 
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From (KTET\f . (KTE^ we get 

In 1/(0)1 > -€o{h)h~\ In \f{w)\ < eo(/i)/i"". (8.185) 

By fl8.18ip . we may assume that the last estimate holds in a larger disc, say 
i?(-(f |^),2ro). Let wi,...,wm be the zeros of / in D{-{f\^),3ro/2). 
Then it is standard to get the factorization 

M 1 

f{w) = e^^^^]J{w-w,), weD{-{^\^)Aro/S), (8.186) 

together with the bounds 

\^g{w)\ < O(eo(/i)/i^"), M = O(eo(/i)/i""). (8.187) 

See for instance Section 5 in [86] where further references are also given. 
For < e < 1, put 

Q{e) = {r G [0,ro[; G -Dq,O q,i such that \w\ = r and |/(u^)| < e}. 

(8.188) 

If r G f2(e) and is a corresponding point in Dq,o „!, we have with rj = \w 



M M 



Y[\r-rj\ < Y[\w-Wj\ < eexp(0(eo(/i)/i"")). (8.189) 
1 1 

Then at least one of the factors |r— rj| is bounded by (ee*^'-'^''^'^-"^ ")'ji/m_ 
particular, the Lebesgue measure A(f2(e)) of f2(e) is bounded by 2M(ee'^'^''°'^'')''~"^)^/*'^. 
Noticing that the last bound increases with M when the last member of 
(183891) is < 1, we get 

Proposition 8.30 Let G C''^ with \a^\ = R and assume that e > is 
small enough so that the last member of l[8.189\} is < 1. Then 

A({rG[0,ro]; \a° + ra^\<R, |F(a° + m^)| < e}) < (8.190) 

Here and in the following, the symbol 0{1) in a denominator indicates a 
bounded positive quantity. 

Typically, we can choose e = exp — 0^ for some small a > and then 
the upper bound in fl8.190p becomes 

exp(C'(P 



^ C(l)/i"- 
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Now we equip Bqd[0, R) with a probability measure of the form 

P{da) = C{h)e'^^°'^L{da), (8.191) 

where L{da) is the Lebesgue measure, $ is a function which depends on 
h and satisfies 

|V$| = 0{h-^^), (8.192) 

and C{h) is the appropriate normahzation constant. 

Writing a = + Rra^, < r < ro(a^), G S'^^~^, ^ < tq < I, we get 

P{da) = C{h)e^^'^^r^^-^drS{da^), (8.193) 

where 0(r) = (j)aO,a^{r) = $(a° + rRa^) so that (j)'{r) = 0{h-^^), A^s = 
N4 + M. Here S'((ia^) denotes the Lebesgue measure on S"^^"^. 
For a fixed a^, we consider the normahzed measure 

^(dr) = C(/i)e'^("V2^-Mr (8.194) 

on [0, ro(a^)] and we want to show an estimate similar to (18.1901) for fi instead 
of A. Write e'^'-^V^'^"^ = exp((/)(r) + {2D — 1) Inr) and consider the derivative 
of the exponent, 

^// N 2D - 1 
r 

This derivative is > for r < ^' ''^ =: 2ro, where we may assume that 
2?^o < ''^o- Introduce the measure Jl > fi hj 

Jl(dr) = C(/i)e'^(^-^''V^^-Mr, r^^^x := max(r,?o)- (8.195) 

Since /I([0,fo]) < /i([ro, 2ro]), we get 

/I([0,r(ai)])<O(l). (8.196) 

We can write 

Jl{dr) = d{h)e'^^'Ur, (8.197) 

where 

^'(r) = C»(l)(/i-^5 + /i-^3+[A^5-iV3]+) = c»(/i-^6)^ (8.198) 
Ne = max{Ns,N5). 

Cf dHUHD. 
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We now decompose [0, ro(«^)] into x /i intervals of length x h'^^. If / 
is such an interval, we see that 

\{dr) Jljdr) \{dr) 

cm -W)- W ^ ^ ^ 

From f l8190|) . IKT99^ we get when the right hand side of (KT89^ is < 1, 
|F(a» + ,./^a-)| <.})/«/) < |ilff) exp(^^l„.) 



n 



- 0(l)'r-^^exp,^.,.) 

Multiplying with /I(J) and summing the estimates over I we get 

/I({r G [0,r{a% \F{a'+rRa')\ < e}) < 0(1)/,-^^!^ exp( ^^^^^"^^^^ Ine). 

(8.200) 

Since fi < Jl, we get the same estimate with Jl replaced by /i. Then from 
(183931) we get 

Proposition 8.31 Lei e > 6e small enough for the right hand side of 
^S.lHfA) to be < 1. Then 

P(|F(a)| < e) < ^(l)^"''^^exp(^^^^lne). (8.201) 

Remark 8.32 In the case when R has real coefficients, we may assume that 
the eigenfunctions ej are real, and from the observation after Proposition 
18.241 we see that we can choose ao above to be real. The discussion above 
can then be restricted to the case of real and hence to real a. We can then 
introduce the probability measure P as in fl8.19ip on the real ball B^d{0, R). 
The subsequent discussion goes through without any changes, and we still 
have the conclusion of Proposition I8.31[ 

8.10 End of the proof of the main result 

We now work under the assumptions of Theorem 18.11 For 2; in a fixed neigh- 
borhood of r, we rephrase fl8.176p as 

|detP5,,| <exp^(-^ J J \n\p,\dxd^ + eo{h)), (8.202) 
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where eo{h) is given in fl8.16p . Moreover, Proposition 18.311 shows that with 
probabihty 

> 1 - C(l)/i-^''-%(/i)e"^'(W) (8.203) 

we have 

I detP^,,! > eexp(-^(-^) //in IpM^d^, (8-204) 



provided that e > is small enough so that 

The right hand side of (Kim) is < 1, Va^ G S^^-\ (8.205) 

Write e = e"^/''", ?= /I'^lnf Then (E?2m\> holds if 

I>Ceo{h), (8.206) 

for some large constant C. ( 18.2031) . 08. 2041) can be rephrased by saying that 
with probability 

> 1 - 0(l)/i^^«-%(/i)e"*^, (8.207) 

we have 

I detP^,,! > exp ^ij^ JJ\n \p,\dxd^-Z}. (8.208) 

This is of interest for ? in the range 

eo(/i)<?'<l. (8.209) 

Now, let r <s f2 be connected with smooth boundary. Recall that < 
K < 1 and that 



(18. 9 p holds uniformly for all z in some neighborhood of dT. (8.210) 
Then the function 

0(^)= // l^b^M^^^ (8-211) 

is continuous and subharmonic in a neighborhood of dT. We shall apply 
Theorem 16.21 with < r ^ 1 constant, to the holomorphic function 

u{z) = det Ps^z- 

Then, according to fl8.207p . fl8.208p we know that with probability 

C(l)eo(/i) ? 
> 1 -TTF^e 8.212 
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we have 

h''\n\u(Sj)\>(j){zj)-e,j = l,...,N, iV x -. (8.213) 
In a full neighborhood of dT we also have 

/i"ln|M(2)| < 0(2) +C7. (8.214) 



We conclude from Theorem 16.21 that with probability bounded from below 
as in fl8.212p we have for every M > 0: 

mu-\0) n r) - 1 A<l>Lidz)\ < (8.215) 
0{1) fe 



r 



+ fi{dT + D{0,r)) 



where /i denotes the measure A(j)L{dz). 

According to Section 10 in [33], the measure ^A(j)L{dz) is the push for- 
ward under p of (27r)~" times the symplectic volume element, and we can 
replace ^A(j)L{dz) by this push forward in fl8.215p . Moreover u~^{0) is the 
set of eigenvalues of Ps so we can rephrase fl8.215p as 

ma{Ps) n r) - ^^vol {p^\m < (8.216) 



vol {p-\dT + D{0,r))) . 



OiX) 

This concludes the proof of Theorem 18. Ij with P replaced by the slightly 
more general operator Pq. 



9 Almost sure Weyl asymptotics of large eigen- 
values 

9.1 Introduction 

W. Bordeaux Montrieux P] has studied elliptic systems of differential oper- 
ators on with random perturbations of the coefficients, and under some 
additional assumptions, he showed that the large eigenvalues obey the Weyl 
law almost surely. His analysis was based on a reduction to the semi-classical 
case (using essentially the Borel-Cantelli lemma), where he could use and ex- 
tend the methods of Hager [32] • 

The purpose of this section is to describe the work of Bordeaux Mon- 
trieux and the author [10] on the almost sure Weyl asymptotics of the large 
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eigenvalues of elliptic operators on compact manifolds. For simplicity, we 
treat only the scalar case and the random perturbation is a potential. 

Let X be a smooth compact manifold of dimension n. Let be an 
elliptic differential operator on X of order m > 2 with smooth coefficients 
and with principal symbol p{x,^). In local coordinates we get, using standard 
multi-index notation, 

|o|<m |a|=m 

Recall that the ellipticity of P° means that p{x, C.) ^ for 7^ 0. We assume 
that 

p{T*X) ^ C. (9.2) 

Fix a strictly positive smooth density of integration dx on X, so that the 
norm || ■ || and inner product (-j ■ ■) are unambiguously defined. Let 
r : L'^{X) — )■ L'^{X) be the antilinear operator of complex conjugation, given 
by Tu = u. We need the symmetry assumption 

P* = TPT, (or equivalents, = P) (9.3) 

where P* is the formal complex adjoint of P. As in [89] we observe that the 
property (19.31) implies that 

p(x,-0 =p(x,0, (9.4) 

and conversely, if (19.41) holds, then the operator |(P + FPF) has the same 
principal symbol p and satisfies (19. 3p . 

Let R be an elliptic differential operator on X with smooth coefficients, 
which is self-adjoint and strictly positive. Let eo,ei,... be an orthonormal 
basis of eigenfunctions of R so that 

Re, = (/i°)2e„ < < /i? < /i° < ... (9.5) 

Our randomly perturbed operator is 

P° = P + g°(x), (9.6) 

where u is the random parameter and 

00 
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Here we assume that a^{u}) are independent complex Gaussian random vari- 

value 0: 

a°~Ar(0,a|), (9.8) 



ables of variance aj and mean value 0: 



where 

(/i?)-^e-(^")^ < a, < (/.?)-^ (9.9) 

3/2 - i 1 , , 

M = , < /3 < -, p > n, (9.10) 

s ^ e I 

where s, p, e are fixed constants such that 

n n n 
— < s < n , 0<e<s . 

Let H'^(X) be the standard Sobolev space of order s. As will follow from 
considerations below, we have g° G H^{X) almost surely since s < p — f . 
Hence G L°° almost surely, implying that has purely discrete spectrum. 

Consider the function F{uj) = argp(cj) on S*X. For given 6'o G S*^ ~ 
R/(27rZ), iVo G N := N \ {0}, we introduce the property: 

No 

P{9o,No): J]|V^F(w)| 7^0on{wG^*X; F(w) = 0o}- (9.11) 

1 

Notice that if P{6q,Nq) holds, then P{6,Nq) holds for all 6 in some neigh- 
borhood of 6q. 

We can now state our main result. 

Theorem 9.1 (JT^) Assume that m > 2. Let < 9i < 02 < 2tt and 
assume that P{6i,No) and P{62,No) hold for some Nq G N. Let g G 
C°°([^i,^2];]0,oo[) and put 

n,,er,o,x = {re''; 9i<9<92, 0<r< \g{9)}. 

Then for every 6 g]0, | — /3[ there exists C > such that almost surely: 
3C{u) < oo such that for all A G [1, C)o[.- 

n r^,,,,o,A) - j^.^oip-\ri,^.^,^,)\ (9.12) 

Here cr(P°) denotes the spectrum and #(^) denotes the number of elements 
in the set A. In 119. the eigenvalues are counted with their algebraic mul- 
tiplicity. 
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The proof actually allows to have almost surely a simultaneous conclusion 
for a whole family of 6*1, 6*2, (?: 

Theorem 9.2 Assume that m > 2. Let Q be a compact subset of [0,27r]. 
Let A'^o G N and assume that P{6,No) holds uniformly for 9^0. Let Q be 
a subset of {{§,61,62)] 6j G 6,^1 < 62, g G ^2]; ]0, cx)[)} with the 

property that g and 1/g are uniformly bounded in C°°{[6i,62]',]0,oo[) when 
{g^ ^2) varies in Q . Then for every 5 g]0, there exists C > such that 
almost surely: 3C{uj) < 00 such that for all AG [1, oo[ and all {g, 61, 62) G Q , 
we have the estimate hO.l^) . 

The condition fl9.9p allows us to choose aj decaying faster than any neg- 
ative power of /i^. Then from the discussion below, it will follow that 
is almost surely a smooth function. A rough and somewhat intuitive inter- 
pretation of Theorem 19.21 is then that for almost every elliptic operator of 
order > 2 with smooth coefficients on a compact manifold which satisfies the 
conditions (19. 2p . (19. 3p . the large eigenvalues distribute according to Weyl's 
law in sectors with limiting directions that satisfy a weak non-degeneracy 
condition. 

9.2 Volume considerations 

In the next subsection we shall perform a reduction to a semi-classical situa- 
tion and work with K^Pq which has the semi-classical principal symbol p in 
dnH). Again, 

K(t) = vol {p G T*X; \p{p) - <t}, t> 0. (9.13) 

Proposition 9.3 For any compact set K G C = C \ {0}, we have 

V,{t) = 0(r), uniformly for z e K, < t < 1, (9.14) 

with K, = 1/2. 

This follows from: 

Proposition 9.4 Let 7 be the curve {re*^ G C; r = g{6), 6 G S^}, where 
0<ge C\S^). Then 

vol(p-i(7 + D(0,t))) = C(t), t^O. 
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This follows from the fact that the radial derivative of p is 7^ 0. 
Using fl9.1ip . we can prove: 

Proposition 9.5 Let 6q G S^, Nq ^ N and assume that P{6o,Nq) holds. 
Then if < ri < r2 and 7 is the radial segment [ri, r2]e*^'^, we have 

vol{p~\-f + D{0,t))) = t 0. 

Now, let < ^1 < ^2 < 27r, ^ e C°°i[Oi, ^2]; ]0, oo[) and put 

^Lo2;r,,r, = {re''; e,<e< 02, ng{0) < r < r^giO)}, (9.15) 

for < ri < r2 < cx). If < ri < r2 < +00 and P(6'j, A'o) hold for j = 1, 2, 
then the last two propositions imply that 

yoip-\drl,^.^^^^^^ + D{o,t)) = o(tVA^o)^ t 0. (9.16) 
9.3 Semiclassical reduction 

We are interested in the distribution of large eigenvalues ( oi P^, so we make 
a standard reduction to a semi-classical problem by letting < h satisfy 

C=^, h^lCr'/^, (9.17) 

and write 

h-'iP' -0 = h^P'. -z=:P + h^ql - (9.18) 

where 

P = /i™pO = ^ a„(a;;/i)(/iD)". (9.19) 

|a|<m 

Here 

a^(x;/i) = C(/i™-l"l) in (9.20) 
aa{x] h) = a° (x) when |a| = m. 

So P is a standard semi-classical differential operator with semi-classical 
principal symbol p{x,^). 

Our strategy will be to decompose the random perturbation 

where the two terms are independent, and with probability very close to 1, 
6Qui will be a semi-classical random perturbation as in Section [8] while 

Wk^H' < h, (9.21) 
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and 

s^]l^.P-l^[ (9-22) 

is fixed. Then /;,™P° will be viewed as a random perturbation of /i^P^+Zc^^. In 
order to achieve this without extra assumptions on the order m, we will also 
have to represent some of our eigenvalues ttj(w) as sums of two independent 
Gaussian random variables. 
We start by examining when 

\\h^ql\\Hs<h. (9.23) 



Proposition 9.6 There is a constant C > such that 119. 23\) holds with 
probability 

Here is a brief outline of the proof. We have 

oo 

/^'"gO = J2 (^j{uj)ej, aj = h'^a^^ ~ Ar(0, {h'^ajf), (9.24) 



and the aj are independent. Now, by the functional characterization of H'^ 
in Subsection 18.21 we get 

oo 

\\h"'^^rHs^Y.\^^^'y<^)\"^ (9-25) 


where {jJ^Yaj ~ A/'(0, (5j)^) are independent random variables and aj = 

Combining this with Proposition 17.51 and standard Weyl asymptotics for 
R leads to the result. 
Write 

ql = i. + ql (9.26) 

= E ll= H (9.27) 

From Proposition 19.61 and its proof, we know that 

Wh'^q'jHs < h with probability > 1 - exp(Co - ^^J^(^)- (9-28) 
We write 

P + h"^q'^ = iP + h'^ql) + h"^ql 
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Theorem 18. II can be apphed with P replaced by the perturbation P + h'^q'^, 
provided that H/i^g^llH^ < h. 

The next question is then wether /I'^g^ can be written as Toh'^'^^'^^q^ where 
qco = J2o<hnO<L^j^j |q^|co < R with probability close to 1. This turns 
out to be impossible without extra assumptions. 

In order to avoid such an extra assumption, we shall now represent a'j 
for hfi^ < L as the sum of two independent Gaussian random variables. Let 
Jo = jo{h) be the largest j for which hfi^ < L. Put 

a' = ^h^e-^^^\ where K > p{M + 1), C> 1 (9.29) 

so that a' < \aj for 1 < j < jo{h). The factor is needed only when 
13 = 0. 

For j < Jq, we may assume that a^{uj) = a'j{u) + a'jiuj), where a'j ~ 
M{0, (o"')^), a'j ~ A/'(0, (o"")^) are independent random variables and 



2 (l\2 I (ll\2 



so that 

< = - i^r X a,. 

Put gi = + g", where 

hfi°<L hii°<L 

Now (cf dnSn])) we write 

P + /.-gO = (P + /i™(g: + g^)) + h"^q'^. 



Theorem 18.11 is valid for random perturbations of 

P,:=P + h"^{q':, + ql), 

provided that \\h"^{q'^ + q'^)\\H^ < h, which again holds with a probability 
as in (19.281) . The new random perturbation is now h^q'^ which we write as 
ro/i^^^^"gaj, where g^ takes the form 

U^)= Yl f^ji^h^ (9-30) 

0<V»<L 

with new independent random variables 

A- = -/I'^-^M^^a' (u;) ~ Ar(0, (l/i-"2^i-V(/i))2). (9.31) 

To To 
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Now, by Proposition 17. 5[ 



To C(l)(/i"^-2^i-"a'(/i))2'' 

Here by Weyl's law for the distribution of eigenvalues of elliptic self-adjoint 
differential operators, we have D x (L/h)"'. Moreover, L,R behave like 
certain powers of h. 

• In the case when /3 = 0, we choose Tq = h^/'^. Then for any a > we 
get 

P(|/3|ci. >i?) <Cexp(- ^ ~ 



for any given fixed a, provided we choose K large enough in fl9.29p . 

• In the case /3 > we get the same conclusion with tq = h~^a' if K is 
large enough. 

In both cases, we see that the independent random variables /3j in fl9.30p . 
fl9.3ip have a joint probability density C{h)e'^^°'''^^ L{da), satisfying f lS.lSp for 
some depending on K. 

With K = 1/2, we put 

e,{h) = h^{{\n]-Y + \n-), 

h To 

where Tq is chosen as above. Notice that €o{h) is of the order of magnitude 
h'^~^ up to a power of In Then Theorem 18.11 gives: 

Proposition 9.7 There exists a constant A'4 > depending on p,n,m such 
that the following holds: Let F d C have piecewise smooth boundary. Then 
3C > such that for 0<r<l/C, ?> Ceo{h), we have with probability 

> 1 Ceo{h) _ ^^-^ /g_32) 

,^^n+max(n{M+l),N4+M) ^ ' 

that 

mh-P^) n r) - ^^vol ip^\m < ^C- + vol ip~\dT + D{0, r)))). 

(9.33) 

As already noted, this gives Weyl asymptotics provided that 

volp"^((9r + D(0,r)) = 0(r°), (9.34) 
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for some a g]0, 1] (which would automatically be the case if k had been larger 
than 1/2 instead of being equal to 1/2), and we can then choose r = e^/(^+"), 
SO that the right hand side of fl9.33p becomes < Cei+"/i~". 

As in [SHI |89] we also observe that if F belongs to a family Q of domains 
satisfying the assumptions of the Proposition uniformly, then with probabil- 
ity 

^ ^ Ct^ _e-^ - Ce"A, (9.35) 

^2/^n+max(n(A/+l),Ar4+A/) 

the estimate fl9.33p holds uniformly and simultaneously for all F G ^. 
9.4 End of the proof 

Let 61,62, Nq be as in Theorem WJ\ so that P{6i,No) and P{62,No) hold. 
Combining the propositions 19.31 19. 4[ 19. 5[ we see that fl9.34p holds with a = 
I/Nq when F = F^^ g^.^ A > fixed, and Proposition 19.71 gives: 



Proposition 9.8 With the parameters as in Proposition 9/7 and for every 
a g]0, -^[, we have with probability 

> 1 ^ — —e C'ow - Ce~^ (9.36) 

J^n+max{n{M+l),Ni+M) 

that 



ma{h-P^) n Tl,^.^,^,) - j^yol ip-\Tl,^,,,,))\ < (9.37) 

Moreover, the conclusion 1^9.31^ is valid simultaneously for all A G [1,2] and 
all {61, 62) in a set where P{6i, Nq), P{62, Nq) hold uniformly, with probability 

>1 Ceo{h) ^g^^g^ 

•^Tn% /^n+max(n(M+l),Af4+M) 

For < 5 < 1, choose e = h~^eo < Ch^-^-^{\nlf, so that e/eo = h~^. 
Then for some we have for every a &]0,1/Nq[ that 

m^{h-P.)nTl,^.^,^,)-j^yo\{^^^^^ 

(9.39) 

simultaneously for 1 < A < 2 and all {61,62) in a set where P{6i,No), 
P{62,Nq) hold uniformly, with probability 

>l~^e-^. (9.40) 
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The upper bound in (I9.39p can be replaced by 

^^(i-/3-25)/(7Vo+l) 

Assuming P{6i, Nq), P{62, Nq), we want to count the number of eigenval- 
ues of Pi^ in 

r — r3 

1,A — J- 6li,02;l,A 

when A — t- oo. Let k{X) be he largest integer k for which 2^ < A and 
decompose 

A:(A)-1 

Tl.A = ( [J r2fe,2'=+i) '-'r2fc(A) ;)^. 


In order to count the eigenvalues of P° in r2fc_2'=+i we define h by /i"*2'^ = 1, 
h = 2"'=/'", so that 

#(a(p°)nr2.,2^+i) = #(a(/i'"P°)nri,2), 

V0l(p"^(r2fc,2fc+l)) = — — — -voi(p~^(ri,2)). 



(27r)" ' ^ (27r/i)" 

Thus, with probability > 1 - C2~^e~'^"' we have 

|#(a(p°) nr2.,2^+0 - ^voip-^(r2.,2^-+0l < ^2^2-^^^-'^-^^)^. 

(9.41) 

Similarly, with probability > 1 - C2^3kix)/^^^-2^>'W/'^/c ^ ^iav^ 

(9.42) 

simultaneously for all A G [A, 2A[. 

Now, we proceed as in [5], using essentially the Borel-Cantelli lemma. 
Use that 

^2^^^e-2'^/^ = 0(l)2^^^e-2'^/^, 

e 

^ 2^-2~™^^"^"^''^W+T = C'(i)A-"-(^"^"^'^)W+T, 

2'=<A 

to conclude that with probability > 1 — C2^^ have 

|#(a(p°) nr2.,,) - ^voip-i(r2.,,)| < CsX—^'^-^-''^^ 
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for all A > 2^. This statement implies Theorem 19.11 □ 
Proof of Theorem 19.21 This is just a minor modification of the proof of 
Theorem 19.11 Indeed, we already used the second part of Proposition 19.71 
to get (19.421) with the probability indicated there. In that estimate we are 
free to vary {g, 61,62) in Q and the same holds for the estimate (19.411) . With 
these modifications, the same proof gives Theorem 19. 2[ □ 



10 Some open problems. 

• The distribution of resonances or scattering poles for certain self-adjoint 
operators like the Schrodinger operator is a very intriguing and difficult 
problem where many basic questions remain unanswered. 

Resonances can be viewed as eigenvalues of a non-self-adjoint opera- 
tor, obtained from the original self-adjoint one by changing the Hilbert 
space. If we take a Schrodinger operator and add a random pertur- 
bation to the potential with support in some fixed compact set, we 
may ask wether with probability close to one the resonances obey some 
kind of Weyl asymptotics. In dimension 1 there is a classical result 
of Zworski |lU2j . saying that we do have Weyl asymptotics without 
any random perturbations, but assuming a non-fiatness condition near 
the end points of the convex hull of the support of the potential. In 
higher dimensions there are results by T. Christiansen [HI |15] and 
Christiansen and P. Hislop [TB] saying that in the "generic case" there 
is (roughly) a lower bound on the number of scattering poles in a se- 
quence of large discs which is of the same order of magnitude as would 
be prescribed by a reasonable Weyl law. The methods of Christiansen 
and Hislop use the analysis of several complex variables and it would 
be interesting to understand the relation with the methods developed 
in our lectures. 

• The damped wave equation in its stationary version is an example of 
a non-self-adjoint operator which is close to a self-adjoint one. The 
eigenvalues are confined to a band parallel to the real axis and since 
the work of Marcus-Matseev [BS] we know that the real parts obey the 
Weyl law with a good remainder estimate. As for the distribution of 
imaginary parts many results are known, G. Lebeau [59], Sjostrand [87] , 
N. Anantharaman [5], S. Nonnenmacher and E. Schenk [80], M. Hitrik, 
Sjostrand, S. Vii Ngoc, [16], [15]. 

The following problem seems to be open: Suppose we add some ran- 
domness to the damping term, and that the underlying geodesic fiow 



108 



is not ergodic (assuming that we work on a compact manifold with- 
out boundary for simphcity). Then with probabihty close to 1, do the 
imaginary parts of the eigenvalues distribute according to a Weyl law 
formulated with the help of the time averages of the damping term? 

What about statistical properties of eigenvalues'? Can it be true, for in- 
stance in the simplest one dimensional situations with Gaussian random 
variables in the perturbation as in [S^, that we have Poisson distribu- 
tion of the eigenvalues? What about correlations between eigenvalues? 
Can we have results similar to the ones that are known for the zeros of 
random polynomials? (Cf [SH |H].) A step in this direction may be the 
recent preprint of T. Christiansen and Zworski [17] where the authors 
study the expectation value of the number of eigenvalues in a domain 
for certain one- dimensional pseudo differential operators with doubly 
periodic symbol. 

It would be of great interest to have the sharpest possible bounds on the 
norm of the resolvent. The general theory of non-self-adjoint operators 
only gives very weak upper bounds which can be sharpened near the 
boundary of the range of the symbol (like for instantance in Section H]). 
In the proofs of the various results on Weyl asymptotics for randomly 
perturbed operators it is quite clear that the random perturbation has 
the effect of improving the upper bounds on the resolvent and some- 
times in dimension one we can even get a polynomial upper bound, as 
was observed by Bordeaux Montrieux [5]. It is currently not clear how 
far these improvements go in higher dimensions. Notice that for ran- 
dom matrices there are results showing that we can have a polynomial 
bound in terms of the size of the matrix. See M. Rudelson [79] . 

Sometimes it is natural to have additional symmetries. For instance in 
the case of the Kramers-Fokker-Planck operator one would like to re- 
strict the random perturbations to the class of such operators. Another 
such class is that of PT-symmetric operators where the question of real- 
ity of the spectrum seems to be of great importance. I recently showed 
that for elliptic PT symmetric operators with PT symmetric random 
perturbations, we have Weyl asymptotics with probabality close to 1. 
In particular most PT symmetric operators have most of their eigen- 
values away from the real axis. 

A long term project may also be to apply the theory to non-linear 
problems. 
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